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Abstract 

In this second paper we complete the classical description of an isolated system of "charged 
positive-energy particles, with Grassmann-valued electric charges and mutual Coulomb interaction, 
plus a transverse electro-magnetic field" in the rest-frame instant form of dynamics. 

In particular we show how to determine a collective variable associated with the internal 3- 
center of mass on the instantaneous 3-spaces, to be eliminated with the constraints Ku n A ~ 0. 
Here fcu n t) is the Lorentz boost generator in the unfaithful internal realization of the Poincare' 
group and its vanishing is the gauge fixing to the rest-frame conditions V(i n f\ ~ 0. We show how to 
find this collective variable for the following isolated systems: a) charged particles with a Coulomb 
plus Darwin mutual interaction; b) transverse radiation field; c) charged particles with a mutual 
Coulomb interaction plus a transverse electro-magnetic field. 

Then we define the Dixon multipolar expansion for the open particle subsystem. We also define 
the relativistic electric dipole approximation of atomic physics in the rest-frame instant form and 
we find the a possible relativistic generalization of the electric dipole representation. 
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I. INTRODUCTION 



In Ref. [1] we gave the description of the isolated system "N charged scalar positive- 
energy particles with Coulomb mutual interaction plus the electro-magnetic field in the 
radiation gauge" in the rest-frame instant form of dynamics. Grassmann-valued electric 
charges were used to regularize the Coulomb self-energies. In that paper there was the 
determination of the regularized Lienard-Wiechert electro-magnetic fields in the absence of 
an incoming radiation field and of the complete potential acting among Coulomb-dressed 
charged particles. This effective potential turned out to be the Coulomb potential plus the 
full relativistic expression of the Darwin potential, which at the lowest order in 1/c 2 reduces 
at the known form of the Darwin potential, till now obtainable only from the Bethe-Salpeter 
approach to relativistic bound states in QFT. This full Darwin plus Coulomb potential can 
be regarded as the classical analogue of the complete transverse as well longitudinal effects 
of the single photon exchange. 

In Ref. [2], quoted as I with its formulas denoted by (1.2.5), we extended the approach of 
Ref. [1]. In particular in I we relaxed the assumption of no ingoing radiation field and we were 
able to show that, at least at the classical level, there exists a canonical transformation from 
an initial situation with a free IN transverse radiation field and decoupled particles with 
a Coulomb plus Darwin interaction, into an interpolating non-radiation transverse electro- 
magnetic field interacting with charged particles having only a Coulomb mutual interaction. 
At later times we can again canonically transform the system to a free OUT radiation field 
plus charged particles with Coulomb plus Darwin mutual interaction. 

These results can be obtained within the formalism of the rest-frame instant form of 
dynamics starting from the approach based on parametrized Minkowski theories of rela- 
tivistic mechanics [3]. The first part of paper I contains a complete updated review of this 
formalism with a resolution of all the ambiguities connected with the relativistic center-of- 
mass notion. As shown in paper I, in the inertial rest frame centered on the Fokker-Pryce 
covariant non-canonical 4-center of inertia Y^{r) of the isolated system, the instantaneous 
3-spaces are the Wigner space-like hyper-planes orthogonal to the conserved 4-momentum 

= Mch^ = Mc (\A + h 2 ; h) of the isolated system, whose invariant mass is M. On 
the Wigner hyper-plane, described by the embedding z^{t,B) = F M (r) + e^(h)a r 1 the 
particles are identified by (Wigner spin-1) 3-vectors ffi(r) determined by the intersection of 
the particle world-line xf(r) = F M (r) + e^(/i)?7[(r) with the Wigner hyperplane. 

In the rest-frame instant form every isolated system is described as an external decou- 
pled canonical non-covariant 4-center of mass x^{r) = (x°(t); xnw ~ -pz x°{t)) parametrized 
(like also F M (r)) in terms of r and of the 6 non-evolving Jacobi data z, h 2 . This decou- 
pled non-covariant free "particle" carries the internal mass M and the rest spin S of the 
isolated system. In the case of charged particles plus the electro-magnetic field they are 



1 As shown in I, (r, a) are observer-dependent 4-radar coordinates with T = t/c being the Lorentz-scalar 

rest time and e%(h) = (—h r ;8 l r h , hr = ) are 3 space-like 4-vectors forming together with W = c%(h) 

l+yl+h 2 

the columns of the standard Wigner boost sending to its rest-frame form Mc (1; 0). 

2 xnw = z/Mc is the classical Newton- Wigner position of the 3-center of mass and we have x^{t) = 
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functions of the 3-coordinates fji(r), their conjugate momenta /^(r) and of the transverse 
electro- magnetic fields A±(t, a), vfj_(r, a) = E±(r, a). There is an external realization of the 
Poincare' group, whose generators depend upon z, h, M and S. 

Inside the Wigner hyper-plane there is an internal realization of the Poincare' group, 
whose generators V(mt) = Mc, V(i nt ), J (int) = <S, K-(int), are evaluated from the energy- 
momentum tensor of the isolated system obtained from the action of the parametrized 
Minkowski theory describing the isolated system in arbitrary non-inertial frames. This 
realization is unfaithful because the rest-frame instant form centered on the Fokker-Pryce 
inertial observer implies (see paper I) the following second class constraints 

V( in t) « 0, iC {int ) « 0. (1.1) 

While the first 3 constraints are the rest-frame conditions (each Wigner hyper-plane is 
a rest frame), the other 3 constraints eliminate the internal 3-center of mass avoiding a 
double counting of the center of mass. Therefore Eqs.(l.l) imply that the dynamics inside 
the Wigner hyper-planes depends only on relative variables satisfying Hamilton equations 
having Mc as Hamiltonian. All the interaction potentials appear only in Mc and JC^ nt ) and 

not in J{i n t) = S, since we are in an instant form of dynamics. 

This complete understanding of the relativistic kinematics of the rest-frame instant form 
of relativistic mechanics was obtained only recently in Ref.[4] and in paper I. The impor- 
tant point is to have a Lagrangian description of the isolated system. This will allow the 
determination of its energy-momentum tensor. Only in this way can we find the interaction- 
dependent internal generators Mc and K,^ n t). In most of the other approaches to relativistic 
mechanics the Lagrangian formulation is unknown, so that we can find only Mc but not 
K,{int)- Therefore we lack the interaction-dependent equations which eliminate the internal 
3-center mass, the gauge variable conjugate to the rest-frame conditions and we are able to 
reconstruct neither the orbits nor the particle world-lines in the rest frame. 

In Ref. [4], even in absence of a Lagrangian, we were able to guess the form of JC^ nt ), 
associated with a given class of invariant masses for a two-body system, so that the internal 
generators satisfy the Poincare' algebra in the rest frame. In terms of the 3- vectors ffi(r), 
Ki(r), i — 1, 2, these generators are (p — ffi — t^) 3 

Mc = y / m2 c 2 + ^( r ) + $ (^( r )) + ^ m 2 c 2 + ^ (r) + $(p2( r )), 

V{i n t) = «i(r) + k 2 (t) m 0, 
-, =* 

J(int) = S = ?7i(t) x «x(r) + 772 (t) x k 2 (t), 

= -fh(r) V / mfc 2 + ^(r) + $(p^(r)) - ff 2 (r) y/m* c 2 + «*(r) + $(/?(r)) « 0, 

(1.2) 



3 Since we have {r)l (r), /% s (t)} = Sij 8 r s , we use the following vector notation: ^(r) = ^[(t)J, K{(t) = 
^K ir (r) = -k[(t)^ . 
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where $(p 2 (r)) is an arbitrary potential 4 . 

If 7r(r) = | [«i(t) — k 2 (t)] is the momentum conjugate to p(r), we have 



Mc « A% nt) c = Y / m2 c 2 + 7 f2( r ) + $(p2(r)) + c 2 + 7r 2 (r) + ^(^(r)). (1.3) 



In Ref. [4] , after having solved the equations /C( int ) ~ and eliminated the internal 3-center 
of mass, we obtained the following orbit and world-line reconstruction 




1 





) 




2 



Ml 



(int) 




Y"(t) + - ((-1) 




) 



4(h) P *(t). 



(1.4) 



^)[^W,^(r)] 



The relative variables p(r) and 7? (t) are to be found as solutions of the Hamilton equations 
with Hamiltonian M.(i n t) c. 

In this paper we complete the classical description of the system of charged particles plus 
a transverse electro-magnetic field. 

In Section II we define the internal collective and relative variables for the isolated sys- 
tem of N charged particles interacting through the Coulomb plus Darwin potential of Ref. [1] . 
Then we will perform the elimination of the internal 3-center of mass and the orbit recon- 
struction in the N=2 case. 

In Section III we define the internal collective and relative variables for the free radiation 
field. Here we rely on the results of Refs. [5] and [6], where for the first time there was 
the study of collective and relative variables for the Klein-Gordon field in the massive and 
massless cases respectively. The adaptation of this results to the rest-frame instant form 
is done in Appendix A. This allows one to eliminate the internal 3-center of mass of the 
transverse radiation field. 

Then in Section IV we will use the results of Sections II and III and the canonical 
transformation of paper I to find the collective and relative variables of the isolated system 
"N charged particles with mutual Coulomb interaction plus the transverse electro- magnetic 
field" defined in paper I and to perform the elimination of the internal 3-center of mass also 
in this case. 

In Section V we review the multipolar expansion of the energy-momentum tensor of the 
particle open subsystem and how to get its (non-canonical) pole-dipole approximation, in 
which the subsystem is viewed as an effective world-line (the 4-center of motion) carrying 
the effective spin of the subsystem. This pole-dipole structure must not be confused with the 



4 The internal boost K.u nt ) satisfying {"P*,/C J } = Sij M in the case of the pure Coulomb interaction Mc = 
^m\c? + ii 2 + ^Jm^c 2 + k 2 + 4 ^ 1 ^ 2 1 are not known. Only in the case of Coulomb plus Darwin potential, 
see next Section, is the internal boost known. 
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pole-dipole structure carried by the non-covariant (Newton- Wigner) external center of mass 
of the isolated system "charged particles with mutual Coulomb interaction plus transverse 
electro-magnetic field" . 

While Section VI is devoted to finding the relativistic version of the electric dipole ap- 
proximation of atomic physics [7, 8, 9], in Section VII we face the problem of defining a 
relativistic extension of the electric dipole representation of the isolated system. We study 
various point canonical transformations suggested by standard atomic physics and we find 
that all of them generate singular terms when the electro-magnetic field is considered dy- 
namical and not an external prescribed field: they are induced by the dipole approximation. 
We then propose a relativistic representation in which the singular terms are replaced by 
contact interactions among the particles: the price for this regularization is that in the 
semi-relativistic limit it is different from the standard electric dipole representation. In 
Appendix B there is the identification of the generating functions of these point canonical 
transformations by the introduction of relativistic Lagrangians for the rest-frame instant 
form description of the isolated system. 

In the final Section there are some concluding remarks and an introduction to the subse- 
quent paper III [10]. 
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II. RELATIVISTIC MECHANICS: THE ISOLATED SYSTEM OF N CHARGED 
PARTICLES INTERACTING THROUGH THE COULOMB PLUS DARWIN PO- 
TENTIAL 

Let us consider an isolated system of N charged positive-energy scalar particles either 
free or interacting through the Coulomb + Darwin potential (1-4.5) (see I and Ref.fl]). 

On the Wigner hyper-planes £ r each particle is described by two canonically conjugate 
Wigner spin-1 3-vectors: the 3-position ^(r) and the 3-momentum /%(r), i = 1, ..,7V, re- 
stricted by the rest-frame conditions V^ n t) ~ (vanishing of the internal 3-momentum) and 
by the conditions eliminating the internal 3-center of mass (vanishing of the internal Lorentz 

— * 

boosts /C(j„ t ) 0). 

Therefore we must introduce collective and relative variables on the Wigner hyper-planes 
and eliminate the collective ones. We consider mostly the N = 2 case, but everything is 
well defined for arbitrary N. There are various possibilities and we must find the one which 
allows one to make the elimination explicitly. 

A. Collective and Relative Variables for N Particles 

In previous papers [11] the problem of replacing the 3-coordinates ^(r), /%(r) inside a 
Wigner hyper-plane with internal collective and relative canonical variables in the rest-frame 
instant form was solved in two different ways: 

a) We may introduce naive collective variables ff + — Yll N Vi (independent of the 
particle masses), k + = V(i nt ) = Z~2l " N ~ and then completing them with either naive 
relative variables p a = J2f=i lai Vi, = ^= J2?=i led «i, a = 1, .., N — 1, or with the 
relative variables in the so-called spin bases [11, 12] 5 . This naive canonical basis is obtained 
with a linear canonical transformation point both in the positions and in the momenta. 

b) We may find the canonical basis whose collective variables are the internal 3-center of 
mass q + 6 and k + = Vu n t) = Yll N ~ 0. To these collective variables are then associated 
relative variables p qa , n ga , a — 1, .., N — 1. However this non-linear canonical transformation 
depends upon the interactions present among the particles (through the internal Poincare' 
generators), so that it is known only for free particles [4, 11]: even in this case it is point 
only in the momenta. The advantage of this canonical transformation would be to allow one 
to write the invariant mass in the form Mc = \J M. 2 c 2 + R 2 + ~ Ai 7 , explicitly showing that 
in the rest-frame the internal mass depends only on relative variables. 



5 Both sets can be used to find the expression of the Dixon multipolcs (see Section V) for a two-particle 
open subsystem in terms of canonical c.o.m and relative canonical variables. 

6 Due to the rest-frame condition V^ int ) ~ 0, we have q + « R + s=s y + , where q + is the internal canonical 
3-center of mass (the internal Newton- Wigner position), y + is the internal Fokker-Pryce 3-center of inertia 
and R+ is the internal M0ller 3-center of energy [4]. These global quantities are defined in terms of the 
internal Poincare' generators [11]. 

7 For two free particles we have Mc = ^Jm\c 2 + n\ + m\c 2 + n\ = J M 2 c? + k 2 + rs Mc = 
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Since it is more convenient to use the naive linear canonical transformation we will use 
the following collective and relative variables which, written in terms of the masses of the 
particles, make it easier to evaluate the non-relativistic limit (m = J2f =1 mi) 

N N 

-> \ — v -> -> try \ T -» 

V+ = Vi, K + = V( int ) = «i, 

i=l m i=l 
N N 

p a = VN V* -f ai ffi, 7T a = —= V] ^ ai Kj, a = 1, .., TV - 1, 

«=i v i=i 

JV-1 

i 

Vi = V+ ' 



VJV a-l 

JV-1 

Ki = — K+ + VN S~] >y ai 7T a , (2.1) 



m 

a=l 



(2.2) 



with the following canonicity conditions 8 

N N JV-1 ^ 

^ ^ Tai 0? ^ ^ lai Ibi &abi ^ ^ lai laj $ij ^5 

i=l i=l a=l 

N m 1 N 

Tai lai / J Tafc) lai ^ ai "T7 / J ^aki 

' m JV 
fc=i fc=i 

AT AT JV-l 

> , — r ai = 0, > 7 a j r bi = 5 ab , \ 7 a j r aj = 5^ — -. 

i=l i=l a=l 

For N = 2 we have 7ll = - 7l2 = ^, r u = y/2 T 12 = 

— * 

Therefore in the two-body case, by introducing the notation 7712 = fj+, k± 2 = k + = V^ n t), 
we have the following collective and relative variables 

mi _ m 2 _ _ _ 

V12 = — Vi H V2, P12 =Vi- V2, 

m m 

^ , ^ n rn 2 _ mi, 

«12 = «1 + ^2 ~ 0, 7Ti2 = «l - K 2 , 

m m 



Vi = V12 H P12, ?72 = ?7i2 P12, 

m m 

Ki = K12+7T12, K 2 = K12-7T12. (2.3) 

m m 



Eqs.(2.1) describe a family of canonical transformations, because the 7 a i's depend on i(A7 — l)(iV — 2) 
free independent parameters. 
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We use the notation m = m 1 + m 2 = m x (1 + 2ja) , fi = = m 2 (1 + ^r)" 1 - For 



m 2 >> mi we have r/ 2 « r/ 12 - Jj£ p 12 , ^ « 7f 12 + p i2 . 



mi • 



The collective variable 7712 (t) has to be determined in terms of pi 2 (r) and 7?i 2 (r) by means 



of the gauge fixings /C(j nt ) = — M _R + 0. 



B. Two Charged Particles interacting through the Coulomb plus Darwin Potential 

From Ref.[l] we get the following internal Poincare' generators for N = 2 9 



£ {int) = M c 2 = cV y mf c 2 + K t + —J; — K- + V DARWIN (fj 1 ( y T) - ?7 2 (t); k^t)), 
i=i 47r|?7 1 -?7 2 | 

— » ~ ~ 

V(int) = «1 + «2 ~ 0, 
2 

i7(int) — x 
i=l 

2 

"tfl2 [E 



(int) 



m| c 2 + Ki + 



i=i 



+ 



+ 



Qi Q2 ( ^ ' % 2 ^ia(«i> «2, P12) - 2 ^± S2 («2, P12) 



( 



+ 



«2 • 



2 y m 2 c 2 + «i 

I <% 12 ^12(«l, «2, P12) - 2 Xj_Sl(Kl, P12)] 



2 y m 2 c 2 + k 2 



)]■ 



""^2 /on ^2 mi / ~2 

- Pl2 I V m l C + K l V m 2 C + K 2 + 

m v m v 



+ 



mi /t 2 • 



m 2 «i 



5 %2 /Cl2 ( /?1 ' ^ 2 ' P12) - 2 4ls2(k2, p 12 ) 



2 m y mf c 2 + /% 

I % 2 ^12(«l, «2, P12) - 2 X]_Sl(«l, P12) 



2 m ^/m 2 c 2 + 



Z.2 
2 



])" 



Vdaw/jv = Vd(??i - ff 2 ,Ki,K 2 ) is given in Eq.(6.19) of Ref.[l] and in Eq.(I.4.5). J( int ) is given in 
Eq.(6.40) and JC {mt) in Eq.(6.46) of Ref.[l]. 
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- Yc ® 2 (\/ m l ° 2 + + / m 2 c2 + «2 % 2 ) ^12(«1, «2, Pl 2 ) ~ 

QlQ2 /" rf3(J / 7T±Sl(^- ^fPl2,^l) + ^±S2^ + ^Pl2»«2) \ _ 



2 / j3 



(f (7(7 



Pi 2 O"-— P12 

m 2 ~ ^ . z, m 1 : 

K±S1\0' Pl2) K lJ " 7T±S2(0" H P12, «2j + 

L m m 

+ B S1 (a - — p 12 , Ki) ■ B S 2{cr + — p 12 , k 2 ) 

= f -McR + ^ 0, (2.4) 

where we have used Eqs. (2.3) in the final expression of the internal boost. 

The above forms for the generators were found in Ref.fl] by first eliminating the elec- 
tromagnetic degrees of freedom by forcing them to coincide with the semiclassical phase 
space Lienard-Wiechert solution by means of second class constraints (no ingoing radiation 
field). Then having gone to Dirac brackets with respect to these constraints, we found 
new (Coulomb-dressed) canonical variables ff i ,for the particles leading finally to a re- 
duced phase space containing only particles with mutual instantaneous act ion- at- a- distance 
interactions. 

These new canonical variables fj^ given in Eqs. (5.51) of Ref.[l], have the following 
expression 



Hr) = *<t) - (-r* \ q, e Qj g/Cl2( " l(r) '^' rfl(r) " r/2(r)) , 

/C 12 (r) = -/C 21 (r) = J d s a [A ±S1 ■ n ±S2 - n ±si ■ A ±S2 ](r, a). (2.5) 



The a-dimensional quantity /Ci 2 (r) is given in Eq.(5.35) of Ref. [1] and in Eq. (1-3.5), while 
the Lienard-Wiechert fields A ±S i, n ±S i are given in Eqs. (1-2.50), (1-2.51). 

Therefore JC^ nt ) ~ can be solved to get ff 12 ~ ^ 12 [/?i, k 2 , p 12 ], so that, by taking into 
account the rest-frame conditions Kyi ~ 0, from Eq.(2.4) we get 



Vl ~ ^fe,^] + Pl2, V2 ~ »7l2 [Pl2,7Tl2] p 12 , «! « -«2 « 7T 12 , 

m m 



with P 12 (t) and 7Ti 2 (t) solutions of the Hamilton equations with Mc as Hamiltonian. 
Then the inverse of the canonical transformation (2.5) allows one to get 



(2.6) 
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m = ~Vi{r)-{-) i+1 \QiY,^ 



9Ki2(ki(t),K2(7-),^(t) -%( r )) 



- ft 2 [ft 2 , *«] + H' +1 ^ ft, - (-)• i Oi E Q, ^ 12( " 12(T) a f 12(T> '^ (T>) . 

* M z mm 1 n V" n ^i2(«i( r )»«2(T),7yi(r) -7f 2 (r)) 

Ki(r) = /ti(r) + - > Qj -5 , (2.7) 

with AC 12 being the same function of its tilded arguments as K 12 and its derivatives with 
respect to 7Ti 2 are taken on the first argument. 

Finally, by using Eq. (1-2.4) for the embedding of the Wigner hyper-planes in Minkowski 



space-time (P M = Mch^ = Mc (v 1 + h 2 ; h) is the total 4- momentum of the 2-body system 
in an arbitrary inertial frame; /i M = e^(h) and e^{h) are the columns of the standard Wigner 
boost sending to the rest-form Mc(l; 0)), we can reconstruct the relativistic orbits and 
the 4-momenta {pf = m 2 c 2 ) of the particles 



( _\i\ n ST n ^^i2(7ri2(r),-7fi 2 (r),p 12 (r)) ^ 
[ j 2 Wl dn„ r )' 



-'7Tl2i 



pf(r) » y/m^ + ^ 12 (r) h" + R i+1 ^ 2 (r). (2.8) 
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III. COLLECTIVE AND RELATIVE VARIABLES FOR THE REST-FRAME DE- 
SCRIPTION OF THE RADIATION FIELD IN THE RADIATION GAUGE 



Let us now consider the rest-frame instant form of a radiation field in the radiation 
gauge, which was given in Eqs.(I-2.33)-(I-2.34) of paper I. Given the internal Poincare' 
generators of Eqs. (1-2.35), we must find a canonical basis spanned by collective and relative 
variables allowing one to eliminate explicitly the internal 3-center of mass as required by the 
constraints V rad ~ and lC rad ~ 0. An approach to the determination of such variables has 
been done in Ref . [5] for the massive Klein-Gordon field and in Ref . [6] for massless fields. The 
rest-frame instant form of this approach was given in Ref. [13] for the massive Klein-Gordon 
field. In Appendix A there is review and the adaptation to the rest-frame instant form of 
the results for the massless scalar field. In this Section we will use these methods for the 
transverse radiation field. 



A. The Radiation Field and its Internal Poincare' Generators 

As shown in I, Eqs. (I-2.33)-(I-2.34), we have the following rest-frame representation of 
the radiation field in the radiation gauge on the Wigner hyper-planes 10 



A=l,2 

d - 

n±rad(T, a) = E± rad (r, a) = - — A± rad (T, a) = 

= i I dku{k) J2 [ax(k)e- l[ ^ )T -^ ] -al(k)e i[uJ ^ T -^ ] 

J A=l,2 

— # —* —* 

B r ad(r, a) = d x A ±rad (r, a) = 



10 „A 



a = (a T = T;a r ); k = (k T = \k\ = w(fc);/c r ), k — , with k Wigner spin-1 3-vector and k T Lorentz 
scalar; dk 



2ui(k) (2tt)3 



, n(k) = 2u(k) (2tt) 3 , [dk] = [l~ 2 ]. 
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a x (k) 



I 



d 3 a e x (k) ■ u(k) A ±rad (r, a) - i 7r ±rad (r, a) 



-i k-a 



{A r ±rad (r, a), 7r s ±rad (r, a ± )} = -cP r ± s (a) 5 3 (a - a 1 ), 



{a x (k),a* x/ (k')} 
{a x (k),a x >(k')} 



-i 2 uo{k) c (2tt) 3 5 xx , 5 3 (k - k') d = - i Q(k) c8 xx , 5 3 (k - k'), 
{a* x (k),a* x ,(k')} = 0, 



S ij = £ e\(k) eft*) + —, k- e x (k) = 0, 



A=l,2 



e x (k) ■ e x ,(k) = 6 XX ,, -|- ■ [ei(fc) x e 2 (k)] = 1, 

\k\ - 



(3.1) 



and the following expression for the internal Poincare' generators of the radiation field [ 

T)A ('T>t c I „ Jl/t „. -p \ sju 1 p urs sjrs sjrs p rsu q-u 1 

' rad \' rad <-radl ^ lvl rad^i ' radji O rad 2 ° radi ^ rad c ^ radi 



M rad c 2 = S rad = cV, 



T 

rad 



i / So [* 



^2 , ry2 

±rad ' rad 



A=l,2 



r,a)=^2 dku(k)a* x (k)a x (k), 



V rad = - I d A a 



J d 3 a 7f ±rad x B rad ^ (r, a) = - ^ j dkka* x (k) a x (k) « 0, 

^ " ° A=l,2 J 



Jrad = S rad = - j d 3 a G X {^l_rad X B rad ) (r, <t) = 



9 



- y dka* x (k)k x — a A (/c) + 

A 



+ Y c J2j d ~ k «y (*) - «a(*) «v (*)] • (* x ^) e v (*) - 

AA 
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^rad 



J, 



rad 



— I d 3 a r 

2c 



'- J dka* x {k)u{k) 



_d_ 



^Irad + B rad ( r > °0 = 

a x (k) + 



+ 



2c 



A,A'=1,2 



a x (k) a*, (k) - a* x (k) a x > (k) e x (k) ■ u(k) 



, T de x >{k) 



dk" 



0, 



h 



rad 



iJrad ' k 

\k\ 



I 



dk 



a* 2 (k) ai(k) — a\(k) a 2 (k) 



(3.2) 



where in the last line we defined the helicity . 

As shown in Appendix A based on Ref. [6] 12 , we must have a\(k), da\(k) E L 2 (R 3 , d 3 k) 
for the existence of the previous ten integrals (and of the occupation number N\ = 
J dk a* x (k) a\(k)) as finite quantities. Moreover one can show [6] the existence of the following 

A — * 3 ^ — * 3 

behavior: i) |a A (A:) | — >|£|_ ) . 00 \k\~ : 2~ p with p > 0; ii) |oa(^)| ^\k\^o \k\^ with e > 0. 

By going to the circular basis a±(k) = -j= (ai(k) =F ia>2(k)), e±(fc) = ^ ( e i(^) ± ^^(A;)), 
and by introducing the modulus-phase variables ([/ CT ] = [m/ 5 t -2 ]) 



= a* a (k)a a (k), <f> a (k) = lln a ° : " : 



a <?(k) = V Ia(k)e 



2 * a*M 

{Uk)AAk')} = ms aa ,sHk-k'), 

the helicity and the Poincare generators of the field configuration become 



(3.3) 



h r ad 
^rad 



- / dk 



[i + (k)-i4k) 

= \ Yl J d~kk A al(k)a a (k) = ^ J dkk A I a (k) = 



= [VI, 



'rad 



rad 



= M rad C; Vrad ~ , 



11 Dimensions: [h] = [J] = [mPr 1 ]. 

12 For the free radiation field the Fourier coefficients and the modulus-phase variables are r- independent. 
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£ad = ~\ J2 j d kl a (k) (kX (f>cr(k) + 

[l + (k) - /_(£)] [e + (£) + e_(£)] .(JxA)[ - + 



+ — I dk 

2c 



+ - y dfce + (fc) x eL(jb) |/ + (Jfe) - J_(Jfe) 
Kad = ~\ E/ dkI a (kMk)-^ ( j )a (k) + 

^ y rffc [/+ ( jfe) - /_ (jfe) e+(fc) + e_ (jfe) 



+ 



"$)-^ r ?-(&)- «0. (3.4) 



As shown in Appendix A and Ref. [6], the existence of finite values of these quantities 
requires [6] : 

i) \I a (k) \ — hu^^ |A;|~ 3 ~ <5 with 5 > 0; \<p a (k) \ — ), |jf|_ ( . 00 |&| (this is due to the transformation 
properties under Poincare transformations; (j)' a (k) = (^(A -1 k) + k ■ a + a #(/c, A -1 )); 



ii) \k\- 2+ ' with e > 0; \<f> a (k)\ 



|k|-0 



\k\ a with a > —5. 



Let us replace the canonical variables I a (k), (f) a (k) with the following new canonical basis 

^rad(k), ®rad(k), \ ra d(k) , prad(k) 13 

^rad(k) = J -W' ^rad{k) = ^(*)> 

a=± a=± 

X rad (k) = I + (k) - I4k), p rad {k) = \ [Mk) - <P4k)] , 

{^rad(k),^rad(k')} = {Kad(k) , p rad $ ')} = ^(k)5 3 (k - k' ) , 



I4k) = - (^rad(k)+\rad(kf), 



I4k) = i (*rad(k) - Xradikf), 



<f>4k) = ®rad(k) + Prad(k), <f>-(k) = <&rad(k) ~ Prad(k), 

in terms of which the helicity and the Poincare generators become 



(3.5) 



13 Dimensions: [* rod ] = [X rad ] = [I a ] = [ml 5 t 2 \. 
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Kad = - / dk\ rad (k), 



^rad 



1 

c 



■'Tad 



rad 



M rad c; P ra d w 



)■ 



J, 



rad 



9 



rfA; rad (fc) ( k x — =, ) $ rad (/c) - 
ok' 



dk 



Kad(k) (k x p rad (k) 



+ 



+ 
+ 



dk' 

/ dk\ rad (k) \e+(k) +e_(fc)l • f fc x — ^) 

^CJ L -I V OK' 



II 



dk\ rad (k) e + (k) x e_(fc) 



d 



Kad = ~ dk^ rad (k)uj(k) -wT-^rad(k) + - / dk X rad (k)uj(k) — p ra d{k) + 



e_(fc) — e + (fc) 



d 



+ 



dk" 



+ 



-I 

2c J 



d 



dk\ rad (k)\e + (k) +e_(A;)J • u(k)—\e-(k) - e + (k) 



0. 



(3.6) 



One has: 



|k|-0 



fc|- 3+e withe>0; \k\ a with a > 1-e; |A ra(i (fc)| 



i) 

with 5 > 0; |p r ad(*OI ^igi^o with ^ > 



'1*1- 



-ol*l' 



-2-<5 



ii) 



|fc|— >oo 



|£|~ 3 -t with 7 > 0; |$(jfe)| 



\k\— >oo 



|fc| ; |A ra(i (fc)| 



|fc|— >oo 



|£|~ 2 -x with 



X > 0; \prad(k)\ ^|fc Hoo l&l' 3 with /3 > x- 



B. Collective and Relative Variables for the Radiation Field 

We can now use the results of Appendix A to replace the canonical variables ^ r ad{k) and 
^rad(k) with a set of collective and relative canonical variables X^ ad , V£, H ra d(k), K. ra d(k). 

To them must be added the remaining canonical variables \ ra d(k) and p ra d(k) describing 
the difference of the modulus-phase variables with different circular polarization (they are 
already of the type of relative variables). 

By using Eqs.(A6), (A7), (A9) and the function T{k, V r B ad ) of Eq.(A8) with uj{k) = k and 
with the functions F T (k) ) F(k) given in Eqs.(All), the collective and relative variables are 
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X. 



rad 



Israel 



oVaa 



dkk A V rad (k) = / dkk A T{k,Vf ad ), 



Hrad(k) = J dqQ(k,q) ^ rad {q) - u(q) F T (q) J dq ± u(qi) ^ ra d(, qi) + 



+ F(q)q- J rfgigi^ rad (gi) 
K rad (k) = V % <$> rad (k), 



^rad (k), Prad(k), 



AB 



\ radi ' radS 

{H rad (k),K rad (h)} = (2n) 3 2Lu(k) 5 3 (k-h) 

{Kad(k),Prad(kl)} = (2vr) 3 2 Lu(k) 5 3 (k - 



(3.7) 



with all the other Poisson brackets vanishing. The operator T>$ and its Green function are 
defined in Ref.[6] and Appendix A. By construction X A ad is a ^-center of phase. 



The helicity and Poincare generators become [D rs = k r 



dk s 



k s ^, D Tr = u{k) 



h 



rad 



dk X r ad(k), 



' 1 



rad 



■'rad 



rad 



M rad C; Trad ~ 



). 



Urad $rad -^rad X T^rad H~ Sradi 

Srad = — / dkH rad (k) (k x K rad (k) - 

/ dk\ rad (k) (k x p rad (k) + 

c J \ dk J 

+ - c I dk\ rad (k) [e_(fc) x e+(k) + % - (e-(k) + e+(fc)) • k x (el(fc) - ?+(£) 



^rad 



Xrad 'Plod X, 



M rad C - 



rad m rad 



dkH rad (k)D Tr K rad 



dk\ rad (k) D Tr p rad (k) - 



(3.8) 
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The gauge fixing IC rad « to the rest-frame conditions V ra d ~ leads to the following 
determination of the 3-center of phase X rad 



X, 



1 



rad 



M rad c 2 



dkH rad (k) D Tr K rad 



(£) + / 



dk \rad(k) D TT Prad(k) + 



+ - / dk\ rad (k) ( e + (k) + e_(k) 



■ D 7 



(k)-e + (k)]) 



(3.9) 



If A ±rad (T,a) satisfies J dkw lm {k) ^ ra d(k) — F(k,Vr ad ) = for / > 2 (see Ref.[6] and 

Appendix A), then we have the following expression of the fields of the transverse radiation 
field 



<7=± 

r A „A 



A ± rad(r, B) = Jd~kJ2 \^$) \j \ V r B ad ) + a X rad $)) + % H rad (k) 



i[k A (X A ad -a A )+a p rad (k)+J dk! dk 2 K ro(i (A:i) Q(k!,k 2 ) A(fe 2 ,fe)] 



+ C.C. 



^±rad(r,a) = i ( dku{k) [^^(^> r B J+^»rf(fc)) +V % H. rad {k) 

e i[k A (X^ ad -a A )+a p rad (k)+f d kidk2K. Tad (ki)g(k u k 2 )A(k2,k)] + c _ c _j ^ 



5 rad (r )C ?) =2 fdk^kx [ea(k)J^(F(k,V* d )+<j\rad(k)) +V^U rad (k) 



a=± 

J[k A (X A ad -a A )+ap rad (k)+J d fei dfc 2 K rad (fei) G(k!,k 2 ) A(k 2 ,k)] 



+ C.C. 



so that the radiation fields have the following dependence upon r and a 



(3.10) 



A±rad{T, (?) — Aj_ ra(i (r — XJ" ad , CT — X mrf ), 
7T±rad(l", (?) = ^ ^radij ~ X^ ad , CT - X rad ), 
B ra d(T,Cr) = B rad {j~ X^ ad , O X rad ), 



(3.11) 



The configurations of the radiation field admitting the collective variables Xf ad have a 

— * — * — * 

"monopole" structure, carried by X rad andp rad ~ 0, plus the "multipoles" K rad (k), H rad (k), 
Kad(k), p r ad(k) (these last multipoles describe the helicity structure). Eqs.(3.9) expresses 
X rad in terms of M rad and of the higher multipoles. 

As in the Klein-Gordon case we have the extra variables X 7 ad , V 7 ad = M rad c with a 
similar interpretation (see Appendix A). 

By using Eq. (1-2.25) we have the following representation of the potential of the radiation 



field: A? ad (Y"(T) + e?(h)ar) = -e^(h) A^ rad (r, a). 
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IV. N CHARGED PARTICLES PLUS THE ELECTRO-MAGNETIC FIELD 

Let us now consider the problem of identifying the collective and relative variable for the 
system of two 14 positive-energy charged particles with mutual Coulomb interaction plus an 
arbitrary electro-magnetic field in the radiation gauge (see Ref.fl] and I). This will allow to 

— * 

solve the constraints K,^ nt ) ~ and to eliminate the internal 3-center of mass like we made 
in the previous two Sections. 

A. The Internal Poincare' Generators before the Canonical Transformation 

From Eq. (1-2.23) we have the following expression of the internal Poincare' generators in 
the original canonical basis ffi(r), «i(r), A±(t, a), 7t±(t, a) 





2 



+ 



Q1Q2 



4tt I fjiij) - fj 2 {r) 





2 



i=l 



em 





2 




The case of N particles follows the same scheme. 
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= - J d 3 aa r T TT (r,a) 



i=i 



C ^mjc^+ «*(r) 



) 



+ 



+ 



Q1Q2 \ V r i(r)+V2(r) 



An c 
d 3 a 



2 

E 



47T 



\#-m(T)\ \a-ff 2 (r) 



+ 



47T C 



+ KZ 



1 

"2c 



\<7-Vi\7)\ 
d 3 aa r (7?i + 5 2 )(r, a) 



«0, 



(4.1) 



For an electro-magnetic field not of the radiation type the quantities £ em , V em , J em -, IC em 
are not conserved and do not satisfy a Lorentz algebra. 

In the canonical basis A±(r, a), n±(r,a), ffi(r), /%(t), it is not clear which are the col- 
lective variables to be eliminated by means of the constraints V(i nt ) ~ and JC( int ) ~ 0, 
because we do not know how to define these variables for an electro-magnetic field not of 
the radiation type. To treat this type of field we now use results obtained in paper I. 



B. The Internal Poincare' Generators after the Canonical Transformation 

After the canonical transformation defined in I, the system is described by a transverse 
radiation field A± rad (r, a), ^± ra d( T ,^), and by Coulomb-dressed particle variables ff^r), 
ki(r), see Eqs. (1-3.6) and (1-3.10). 

From Eqs.(I-4.4), (1-4.1), (1-4.2), (1-4.6) [and (1-3.9) for % and t i3 \ we have that each 
internal Poincare' generator becomes the direct sum of the radiation field one of Eq.(3.2) 
plus the particle one of Eq.(2.4) [see Eq.(6.46) of Ref.[l]] 



i=i 4?r l Vi(V -r) 2 {T)\ 

+ VDARwiNiKiir),^^),^^) -ff 2 (r)) + 



(^matter '^'radj C i 
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2 ^ . 

V(i n t) = ^2 ^( T ) + ~ / d3(T (^±rad X B rad j (r, a) 
i=l ° 



^ 'matter T^rad ~ 0, 



J(mt) = Y1^ X \ j d3(J ° X (*±rad X ^rad) (t, 5) = 

— * — * 
\J matter \Jradi 



1 



1 Q,Q 2 r, -ft 

Vi 1 1- 

1,2 



2 c 



m 2 c 2 + k 1 



- 2 ; U 2a ± 5i(«i,pi 2 jj 



m| c 2 + k 2 



- 1 ^ ^ 2 m 2 c 2 + ^ JL + ^ m 2 c 2 + £ 2 \ ^ ^ _ 
* C \ O Ki a K 2 / 

Q1Q2 f f nxsiia-fj^K!) jr ±S2 (a -ff 2 ,K2) \ _ 

- Q 1 ® 2 J d 3 aa (n± S i(a - ff l7 k ± ) ■ TT± S2 (a - f} 2 , k 2 ) + 
+ B S i(& -f]i,Ki) • -B 52 (<t - ?7 2 , k 2 )) - 

^matter ~\~ K^rad 



£{int) R+ — — - (^matter + ^rad^j R+ ~ 0. (4.2) 



def _1- - I 

C 



Let us remark that in /C( int ) ~ of Eq.(4.2) all the particles terms depend on p 12 (r) 
-%( r ) or a-ffiir). 



21 



We also added the connection of /C( int ) with the M0ller internal 3-center of energy R + 15 . 

Even if all the internal generators are the direct sum of the generators of the two subsys- 
tems, the two subsystems are connected by the rest-frame conditions and by the vanishing of 
the internal boosts, i.e. by the necessity of eliminating the position and momentum of the 
internal overall 3-center of mass. 

C. Semi- Relativistic Expansions of £^ n t') 

and /C( int ) before and after the Canonical 

Transfer mat ion 

The semi-relativistic limit of the generators £(m<) and JC^ nt ) of Eqs.(4.1) is 



c ^-^ rrii c z 8mf 

t=i i=i 1 

c d 2rrii rrii 



fc\int) cJC r Ga[Uei + 0(c *) PS 0, 

AT AT 

/CcaH/ei = - ^ mi $(t) = - ( ^ mi ) ^( n ) ~ ' ( 4 ' 3 ) 
i=l i=l 

where X( n ) is the non-relativistic center of mass, which emerges as the limit of the M0ller 
internal 3-center of energy R + = —dC^ n t)/S^ nt ). See Section G of paper I for the non- 
relativistic limit of the relativistic variables ffi(r), K-iij) and of the relativistic external center 
of mass. 

After the canonical transformation S^ nt ) is given by Eq.(4.2). The semi-relativistic limit 
of its matter part is 



We have R + « ? + « y + due to the rest-frame condition V( int ) ~ 0, where q + is the internal 3-center of 
mass and y + the internal Fokker-Pryce 3-center of inertia. 
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matter 



1 
i=l 



2 ( T ) 



Q1Q2 1 Kir) 

2mi ' 4tt|^ 1 (t) (r) | c 2 8mf 



+ 



+ 



Ql Q 2 



Vl( T )-V2( T ) 
~rf 2 (r)\ 



m 1 m 2 c 2 167r|7h(r) - ff 2 (r) 



-)+0(c- 4 ) 



(4.4) 



By using Eqs. (1-2.51) - (1-2.53) for the Lienard-Wiechert fields we get /Ci 2 = 0(c 3 ), so 
that the semi-relativistic limit of JC^ nt ) of Eq.(4.2) is 



(int) 



= -yjlM( m ,c + A-) + 0( c-) 

i=l 

~ Tel*"* ( j ?lra d + B 2 rad )(T,a) = 



~ * ~ * ~ * 

= ^matter + fcrad = cK Galilei + O(-) « 0, 



^Galilei = -(%+%) X(n) = ~ m i Vi ~ 

i=l 



(in the rest frame). (4.5) 



D. The Collective Variables after the Canonical Transformation 



The results of Section II and III allow us to find the collective and relative variables 
of the matter and transverse radiation field subsystems in the canonical basis ff^r), Ki(r), 

A±rad{T, a), TT±rad(T, ff). 

A) The field quantities A± rad (r, a), n± rad (r, a) are canonically equivalent to the canonical 
basis X T rad , V T rad = M rad c = ^f-, X rad , V rad , Kad(k), p ra d(k), H rad (k), K rad (k), given in 
Eqs. (3. 7), which contains the collective variables for the radiation field. 

B) For the particle subsystem we make the canonical transformation (2.3) on the 
Coulomb-dressed particle variables 



- mi - m 2 u % u u 

V12 = — Vi + — Pl2 = J 7l- J 72> 

m m 

«12 = «1 + «2 ~ 0, 7Ti2 = Kl — «2, 

m m 



% u m 2 - 

»7l = ^12 + Pl2> 

m 

«1 = «12 + 7I"12, 

m 



^2 — ^12 ~ Pl2> 

. m 2 - % 

K 2 = K 12 — 7Ti 2 . 

m 



(4.6) 
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Eqs.(4.2) imply V matt er = «i2 ~ -Prod = - J d 3 a (n ±rad x B rad )(r, a). 

C) Let us now put ff 3 = X rad , ^3 = V rad , with X rad and P ra ^ given by Eqs. (3.9) and 
(3.8) respectively and let us combine it with fj 12 and k 12 of Eqs. (4.6) to get the canonical 
transformation 



Vl2 


Pl2 


% 


«12 


7?12 


«3 



0) 



Pl2 



7T12 



P(12)3 



7r (12)3 



I 



P(12)3 



2 (^12 + ^73) 5 « = «12 + «3 = ^(mi) ~ 0, 



Vl2 ~ V3i 



7T(12)3 



(«12 - K 3 ), 



so that we get 



Vl2 
V3 



V+ 2 P(12)3» «12 = 2 ^ + ^(12)3 ~ 7T(12)3, 

. 1 . . 1 . . . 

V - 2 P(12)3> K 3 = ^ K ~ 7r (12)3 ~ 7T(i 2 )3, 



(4.7) 



Vi 

V2 

Y 

^-rad 

«1 

^2 
Vrad 



i, U . m 2 i. 

^ + 2^(12)3 + —P12, 

- 1 - mi ^ 

^+ 2^(12)3 - ~Pl2> 
V3 — V - 2 ^(12)3) 



m, i 

m 2 
m 2 ,1 



« + 7T(12)3) + 7T12 



mi - - 
7T(12)3 + 7T12, 

m 

m 2 ±, 



I 9 K + 7r (12)3 J - 7T12 ~ 7T(12)3 ~ 7r 12, 

777- ^ 777* 

«3 = ^ K _ 7r ( 12 )3 



-7T 



(12)3- 



(4i 



Therefore our final canonical basis contains the collective variables ff, k and the 
relative variables p 12 , 7fi 2 (depending only upon the particles), K mc ;(/c), H ra d(/c), \ rad (k), 
Prad(k) (depending only upon the "multipoles"the transverse radiation field) and p*( 12 ) 3 , 

^(12)3 (depending on the relative motion of the particle collective variables with respect to 
the radiation field collective variables, i.e. its "monopole" aspect). 



The overall collective variable ff is the natural variable for the solution of the gauge fixings 
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E. The Internal 3-Center of Mass ff from the vanishing of the Internal Boosts after 
the Canonical Transformation. 

We have to find ff of Eq.(4.7) from the vanishing of the internal boosts /C( int ) m in the 
form (4.2), put it into Eqs.(4.8 ) and make the inverse canonical transformation (1-3.10) to 
find the original ffi(r) and the particle world-lines xf (r) like we made in Eqs.(2.8). 

In the boosts in the form (4.2) the first four lines of the particle terms depend on p 12 = 
77 x — 77 2 and 7T12 = ^ K\ — ^ k 2 , but they also get a dependence on 7T(i 2 )3 = \ (ki + K> 2 — V ra d) 
through the dependence on the particle momenta. In the next two particle terms we must 
shift the integration variable to reabsorb the quantity ff + \ f>(\2)z present in a — ff { . 

Finally the last line of the boosts in Eqs.(4.2) is just JC rad of Eqs.(3.8). 

Therefore the internal boosts in the form (4.2) has the following form in the new canonical 
basis (ki(t) are given by Eqs.(4.8)) 



K. 



(int) 



X 



V + g ^(12)3 
■? ( 1 9^12 



2 

E 

i=l 



2 o ^2 1 Qi Q 2 



x 



(ki,k 2 ,p 12 ) 



dPi2 



2 A ± S2{K2, Pn, 



+ 



m\ c 2 + k 1 

i ( 1 9^12(ki,K2,Pi 2 ) 94* ^ \ 

12 ' 1 2 2A ± 5i(«i,Pi 2 J 



+ 



m 2 , c 2 + k 2 



— p 12 — V mi c z + k, V m c z + k 2 + 

V /// m 



+ 



1 Q1Q2 



m 2 
m 



. 1 9ACi2(bi,B2,Pi 2 ) 
2 ^2 



2 A ±S2 (K2,Pl2) 



ml c 2 + /% 
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mi 
m 



. 1 dfci2(Ki,«2,Pi 2 ) o A % \ 

K 2 ■ \ 2 ^ ^^±Sll^l,Pl2j 



m| c 2 + k 2 



])- 



1 Q1Q2 

2 c 

4?T c 

Q1Q2 



~2 (9 / " 2 \ " * * ■» 

m 2 c 2 + + ym^c 2 + K 2 — ^ ) /Ci 2 (ki, « 2 , P12) - 

T ^ !<? + — A2I 



I -? 1 mi r» 

l^ — P12 

C + »7 + ~ P(12)3 



+ 



m Pl2l 



^LSl^ p 12 , «l) • 7T±52(0 J H p 12 , ^2) + 



\ m m 

+ B S1 (CT p 12 , Ki) ■ B S 2{0- H p 12 , K>2) - 

m m J 

~ ^rad ^(12)3 ~ M rad C (ff - - P(i 2 )3^ + "^(H m d, K rad , X rad , p rad ) PS 0, 



where the term D(H ra d, K ra ^, A ra d, p mc z) has the following expression 



(4.9) 



D (H mc f, K ra( ^, \ rad , p r ad 



>--§/ 



rfA; i7 rad (fc) _D or ^(fc) — - j dk \ rad {k) D or p rad (k) 



dk\ rad (k) e+(k) + e-(k) 



■D c 



e-(k) — e+(k) 



(4.10) 



The solution of Eq.(4.9) is 



1 



rj 



A + M rad c 



- - (A - M rad c) p (12)3 - 



m2 / 2 2 , , m l - X2 m l / 2 2 , m2 - \2 , »1 - /? , 
K\m\& + (7Ti2 H 7T(i2)3r \ ^2° + (^12 7T(12)3r + D Pl2 ~ C + 

my m m V m J 

+ ^^(12)3 + (4-11) 

where we have introduced the following definitions («i(r) are given by Eqs.(4.8)) 
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2 

A(p 12 , 7?12, 7T(12) 3 ) = \/ m l c2 + «i + ^ 9l9l 



X 



X 



1=1 



m\ c 2 + k x 



+ 



+ 



. 1 9ACi2(ki,K2,Pi 2 ) o 4 
K 2 ' I 9 Z 71 



±51 



(«l,Pl 2 )) 



Wlj C 2 + K 2 



QlQ 2 



« 7r_L5i((T p 12 , /ti) • 7Tj. S2 (o- 

V m 

+ .651(0- p 12 , ki) • .652(0- H p 12 , Kfcj , 

m m J 



mi 



m 



Pl2»«2) + 



^(Pl2 ^12,^(12)3) 



1 QlQ2 



m 2 
m 



9f>i2 



2A 



±S2 



(«2,Pl2)) 



mi 
m 



«2 



9ATi2(ki,K2,Pi 2 ) 
<9pi2 



-2A 



±Sl(«l 



>Pi 2 )) 



mf c 2 + /% 



m 2 c 2 + k 2 



C(Pl2, 7Ti2, 7T(i2) 3 ) = ( y m 2 C 2 + K 1 — 5- + \J m\ C 2 + K 2 — 5- ) K, 12 {k U « 2 , P 12 ) ~ 



Q1Q2 f d 3 a rf±si(0- gp^jO | ^2(0 + ^Pi 2 ,/t 2 ; 



/ da a 7r_L5i(cr Pi2,«i) -^±52(0 H Pl 2 ,«2j + 

c J \ m m 

+ 5 S l(<7 Pi2,«l) • ^52(0 H Pl2>«2))- (4-12) 

m m / 



Since, as shown in Eqs.(4.2), we have /C(j nt ) = ciCcaiiiei + O(^) ~ we have 77 = 
£(n) + 0(c~ 2 ) ~ in the non-relativistic limit, namely we get the non-relativistic rest frame 

^ Ga iu e i = -mX(n) = ~ J2i=l m i Vi ~ °- 

By using Eqs. (4.11), (3.8), (4.8) and the rest-frame condition Vu n t) = k(t) ~ 0, we 
get the following expressions of the internal invariant mass and of the internal angular 
momentum in Eqs. (4.2) 
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(int) 



Mc 2 w M raa c 2 + 



+ c^Jmfc 2 + (7Ti2 + ^ 7r ( i 2 ) 3 ) 2 + c y m 2 , c 2 + (7fi2 - ^ 7r(i2)3) 2 + 
Q1Q2 . T/ /i i . mi i i m 2 i, s 

+ -5 — " + V D ARWIN\P\2, Kl2 H 7T(12)3; _7r 12 H 7T(12)3j) 

47r|p 12 | m m 



=>2 1 R 2 
^irad -°rad 



(4.13) 



J(mt) = S W p 12 X 7Ti2 + P(i2) 3 X 7T( 12 )3 + S'rod, 
»Srod 



1 /* 1 

-- / dkH rad (k) k x — K rad (k) - 



dk 



d 



+ 



dk X r ad(k) k x — p p rad (fc) + 
^ y dk\ rad (k) e-(k) x e+(fc) + ^ (e_(fc) + e+(fc)) • fc x A (e_(£) - e+(Ai) 



(4.14) 



Eqs.(4.7), (4.8) and (4.9) imply 



^ = t+ ^(12)3 - ^ + M radC ( M «- C ^)3 - 

m2 / 2 2,/- 1 7711 - A2 mi / 22,/- m 2 7 37 . c 

— Wm c 2 + (tt 12 H 7T(i 2 ) 3 ) 2 a m z 2 c 2 + (vr 12 Tr^p) 2 + o 

m V m m V m 



Pl2 ~ 



- C + XL. 7T 



rod ^(12)3 



1 / 



mi / 2 2,/- m 2 7 37 . B 

H < m 2 c 2 + (tt 12 vr (12 ) 3 ) 2 + B 

m V m 



V2 



Mr-„AC V 



.4 + M ma c 



M rad C P(i2)3 — 



/; 12 - C + X r r ad vf (12)3 + D ) , 

m V 



Wmf & + (7Ti2 H 7T(i2) 3 ) 2 + 

m V m m 



- \\m\c 2 + (tTi2 - — 7T(12)3) 2 



B 



k 2 -c + x; a J {12) z + v). 



(4.15) 



These equations, together with the inverse canonical transformation (1-3.10) and by using 
Eqs.(I-2.51), (1-2.52), (1-3.9), (3.8) and (4.8), allow to get the original variables ffi(r) and 
the world-lines (1-2.18) of the particles only in terms of relative variables 
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aftr) = Y»{T) + e?(h)Til 



rjl(r) * ft (r) - ^ ^ 



Qi m d 



C mi <97T(i2)3, 



(i 3 ( 



■ Al_Sl ( s - 7T 



1 ^ 



- A 



m 2 . mi - - \ 

2 P ( 12 ) 3 "^ Pl2 '^ 7r(12)3 + 7ri2 J" 
1 



_Lmd - X^ ad , (7 + - p (12)3 

- / 1 - m 2 - mi - 

•7T±S1 I C - - P(i2) 3 Pl2> 7r (12)3 + ^12 

Qi Q 2 m d f 3 



2 c mi 9 #(12)31 



±51 [<T ~ - p (12)3 - — Pl2, — 7T (12)3 + ^12 J 

/ 1 



- / 1 - m 1 ^ m 2 ^ - \ 

•VT±S2 ^ - 2 ^(12)3 + — Pl2> — 7T(12)3 - ^12 ) ~ 

- U m 2 ^ mi ^ \ 
- vr ±5 i {a - - p (12)3 - — p 12 , — 7r ( i 2 )3 + ^12 J • 



•A 



1 , mi - m 2 - ^ \ 

1 S'2 ( <T ~ ~ P(12)3 + — Pl2. — 7r (12)3 - 7T12 J 



Q 2 m d 



■A 

- A 



c m 2 dTT {12)3 
( - 1 



d 3 a 



n±rad(r - X T rad , p (12)3 ) 



(_ 1 - mi - m 2 - ^ \ 

■±S2 [a - - p (12)3 + — Pl2 , — 7T( 12)3 - 7Ti2 j - 



•71"±S2^0" - ^ P(12)3 + — 

QiQ 2 m d 



m 2 - - \ 

I Pl2> 7r (12)3 - 7T12 + 

m m /J 
/ ,3 [? 1 m 2 i, mi 

/ <l <T A ±S il(J - -P(i2) 3 Pl2, 



VI V2 "* 

2 c m 2 97T(i 2 ) 3j 



7T(12)3 + ^12) 



1 

2 

1 S. 



, , ^ ^ , mi - m 2 - - \ 

■K±S2 [O- - - p (12)3 + — Pl2 , — 7T(i 2 )3 ~ 7T 12 J - 



VT±51 [O- - - P( 12)3 - — Pl2> — T(12)3 + 7T12 J 



V 2 v ; m m 
r h . mi - m 2 

■A±S2 I Cr - - P(i2) 3 H P12, 



7?(12)3 - 7?12j 



(4.16) 



F. The Final Relative Canonical Variables 



The final independent canonical relative variables are 
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a ) Pi2( r )> ^12 (t) (relative motion of the two particles); 

b) P(i2)3( r )! ^(12)3 (t) (relative motion of the "particle system" with respect to the "radi- 
ation field system" ) ; 

c ) X rad( T )> ^ > lad( T ) = \ £rad(r) = M rad (r) c (the energy of the radiation field and its 
conjugate temporal variable, the only surviving collective variables); 

— * — * 

d) H rad (k), ~K rad (k) (relative multipoles of the radiation field with respect to its monopole- 
like collective variables); 

— * —* 

e) Kad(k), Prad(k) (relative variables describing the helicity degrees of freedom of the 
radiation field). 

These variables satisfy Hamilton equations having S( int ) of Eq.(4.13) as Hamiltonian. 
There are the following constants of motion: 

A) The relative 3-momentum 7T(i 2 )3 = \(k>i + k 2 — V rad j ~ —V rad ~ K\ + k 2 : 

{^(12)3, £(int)} = 0. 

The constant of the motion 7T(i 2 )3 connects the particles and the radiation field. Its 
conjugate variable satisfies the kinematical Hamilton equation 



^P(12)3( T ) o_ d£(i n t) _ c m l ^12 + ^(12)3 

dT d ^ ' m a/^ 2 + (^12 + ^(i2) 3 ) 2 



+ c m 2 *12 ~ g ^(12)3 + 

m ^^ + (^-^^3)2 

+ 9 VdARWIn(P 12 ] 7T12 + ^ ^(12)3i ~^12 + g ^(12)3) ^ ^ 

^7T(12)3 



B) The energy of the radiation field: {P^ ad ,£(i n t)} = 0. Its conjugate variable X^ ad (r), 
iKad^rad} = -1, satisfies 

^jf= = (4-18) 



Like for the Klein-Gordon case of Appendix A, the two second class constraints V^ ad — 
E ra d ~ 0, X^ ad — t 0, eliminating the last two collective variables, select a symplectic 
sub-manifold of the surface of constant energy V^ ad = E rad of the radiation field. 

C) The multipoles ii rad (k), K rad (k), X rad (k), p ra d(k). 

Therefore the two subsystems (particles and radiation field), although not coupled in the 
equations of motion due to £(i n t) = ^matter + £ ra d, are nevertheless effectively interacting 
through the rest-frame constraints and their gauge fixings as clear from Eq.(4-17). 
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The canonical variables q + , R + = k 12 , p q , n q , describing the canonical (Coulomb-dressed) 
Newton- Wigner 3-center of mass and the relative motion of the 2-particle subsystem (see 
Subsection IIA), could be found by using Eqs.(4.7) and (4.8). Then we could re-express 
£{int) of Eq.(4.13) and iTjmt) °f Eq.(4.14) in terms of them. In particular these variables 
would allow to write £(i n t) of Eq.(4.13) in the form 



£(int) ~ cy/Mlc 2 + TT {12)3 + ( function of p q ,n q , 7?(i2)3 ••••) +K ad , 



M c = \J mjc 2 + tf q + yj m\c 2 + fi q . (4.19) 

In this way the 2-particle subsystem is visualized as an effective pseudo-particle (an atom 
after quantization) of mass M. plus interactions. However a drawback of these variables 
is that the Coulomb interaction depends upon p 12 and not upon p q , variables which do not 
coincide for k 7^ 0. 



G. Coming back with the Inverse Canonical Transformation 

If we use the inverse canonical transformation (1-3.10) of I, Eqs.(4.7)-(4.8) and the con- 
sequences (4.11)-(4.15) of lC(i n t) ~ 0, we can get the expression of V em , J em and V T em of 
Eqs.(4.1) in terms of the radiation field and of the Coulomb-dressed particles from the com- 
parison of the internal Poincare' generators before and after the canonical transformation. 

Eqs. (4.1) and (4.2) imply 

V ern = J d 3 <7 (lf ± X B) (t, o) « 

s ^ Qi dfj(T) 1Q 1 Q 2 ( d d w 



i=l 



-7T(12)3 + O(-), 

c 



(4.20) 
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Jem = - I d 3 a a x (V ± x Bj (r, a) 

Srad - if) - ^P(12)3) X ^(12)3 ~ 



2 



i=l 



1 gig 2 

2 c 



»7i x 



E 

i=i 



^ hKj X ^ 



dffi 



+ Ki X 



)- 



Srad + O(-), 
C 



(4.21) 



>r _ 1 



em 



2c 



rfV (7ri + J B 2 )(r )( 7) 



Qi Q2 Q1Q2 , l, 7 t u u % \ 1 

I "5 — "J T^— T + - VDarwin{Pl2j ^12, ^(12)3) + 

Anc\p 12 \ 4ttc |pi 2 | c 

Qi Kj ■ A ± (T,ffi(T)) _ 

= ^ / d 3 a[nl rad + B 2 rad )(r,a) + 0(^). (4.22) 



This is in accord with Eq. (1-2.49) whose semi-relativistic limit a,em\{ T -,k) = 
a x (k) e- iuj ^ T +u(k) e x (k)-J2 2 i=1 ^ sli^ll % x (k^t) x fc) +0(c~ 3 ) says the Fourier coef- 
ficients of the transverse electro-magnetic field differ from those of a transverse radiation field 
for particle terms of order 0(c~ l ). This is the error when we replace the electro- magnetic 
field with a radiation field, as is often done in atomic physics. 
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V. THE MULTIPOLAR EXPANSION OF THE PARTICLE ENERGY- 
MOMENTUM TENSOR 



As said in Section I of paper I, in the rest-frame instant form the isolated system of 
charged particles plus the transverse electro-magnetic field can be described as a decoupled 
non-covariant center of mass carrying the invariant mass and the spin of the isolated system, 
i.e. some type of pole-dipole structure. The invariant mass and the spin are evaluated by 
means of the energy- momentum tensor T ab (t, a) = T^ tter {r, a) + T e ^f (r, a) determined by 
using the action (1-2.1) of the parametrized Minkowski theory. 

Let us now consider the open subsystem composed only by the particles, whose non- 
conserved energy-momentum tensor is T^ ter (r, a). The study of its multipolar expansion 
(see Ref.fll]) allows one to replace the extended subsystem with its multipoles when some 
analyticity conditions are satisfied and then to define a pole-dipole approximation of the 
subsystem. As shown in Ref.fll], this can be done also for single strongly bound groups of 
particles, to be used to describe atoms after quantization, inside the particle subsystem by 
identifying the effective energy-momentum of each group. 

Till now this is the only description of a relativistic (either free or interacting with the 
environment) composite system like an atom as a collective point-like system endowed with 
multipolar properties. The lowest approximation of the composite system is the pole-dipole 
approximation, in which the system is simulated with a point-like particle (monopole) of 
mass M c moving along the world-line of an effective 4-center of motion and carrying the 
spin dipole J c (to be used for the magnetic dipole moment). In the rest-frame instant form 
the 4-center of motion has the world-line w£(t) = F M (r) + e^(h) Q(t) and 4-momentum 
PM = WM c c + e>t(h)V r c . 

Let us remark that this pole-dipole approximation is not the "dipole approximation" of 
the semi-relativistic atomic physics, which will be studied in the next Section. 

As shown in Ref.fll], given a non-isolated cluster of particles the main problem is the 
determination of an effective 4-center of motion described by 3-coordinates Q(t) an d with 
world-line w%(t) = F M (r)+e^(/i) Q(t) in the rest frame of the isolated system 16 . The unit 4- 

velocity of this center of motion is u%{t) = w£(t)/\J 1 - C c (r) with w%(t) = h^ + e^ih) Q(r) 
(from paper I we have h» = u"(P) = P^/Mc). By using Sz^(r,a) = e»(h) (a r - C(r)) we 
can define the Dixon multipoles of the cluster with respect to the world-line w%(t) 

e- r " A V) = / d 3 a [a ri - Q(r)lV n ~ C"(r)] T^ tter (r, a). (5.1) 

The mass and momentum monopoles, and the mass, momentum and spin dipoles are, 
respectively 



See Ref. [11] for a discussion of the main choices present in the literature. 
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q T J = M c , q r J = V r c , 

i T = -K - M c C(r) = M c (R r c (r) - C(r)), <£" = f c u {r) - Q{r) V u c , 
P7 = /^^: to (r,a>^[(r)<(r)- 

J i=l 
2 „ 

- Qi d^c{a-^{T)){d r Al + d u Al){T,a), 

i=l J 
2 

P?+PT = E (^W«?(r)+^(r)<(r))- 

i=l 

2 „ 

- 2^Qi d s ac(a-fJ t (T))(d r Al + d u A r ± Yr,a) 

i=i 



ru ur def ruv q- v 

Pc Pc ~ t J c i 



UTT\ 



S? = [e?(h) K" - e v r {h) h?\ q r c TT + e?(h) e^h) {q™ - q, 
= [${h) h v - e»(h) h"] M c (R r c - C) + 
+ e?(h) [e™ T c - (£ V u c - Q V r c ) 

=► r< (P) = K = -e?(h) q r J\ (5.2) 

where J c and JC C are the angular momentum and the boost associated with the chosen 
definition of effective center of motion. 

As shown in Ref.[ll], it is convenient to choose the center of energy R c = —)C C /M c c 

— * — * 

as center of motion. R c and V c (momentum monopole) are the non-canonical variables 
describing the monopole, i.e. the collective pseudo-particle of non-conserved mass M c (mass 
monopole). 

On the world-line of this collective pseudo-particle there is also a spin dipole (the anti- 
symmetric part of the momentum dipole) and the symmetric part of the momentum dipole 
17 . See Ref.fll], Eqs. (5.23) - (5.26). 

Since we are in a Hamiltonian formulation, the constitutive relation between monopole 

3- velocity dC '^ and monopole (non conserved) 3-momentum V c must not be added by hand, 
but can be determined for each choice of collective center of motion [11]. 

As shown in Ref.fll], the total 4-momentum P£ = e^(h) qf T = W M Q + e%(h) V r c and the 
spin dipole tensor S£ v obey the Papapetrou-Dixon-Souriau equations. 



17 This symmetric part of the momentum dipole depends on the electromagnetic potential and is connected 
with the electric dipole: it is directed along the associated relative momentum if one introduces a suitable 
definition of a relative variable. 
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A. The Multipolar Expansion after the Canonical Transformation. 

Let us consider the open subsystem formed by the particles. 
From Eqs.(4.13) and (4.14) we get 

Smatterir) = cJm\c 2 + (^12 + — 7?(i2)3) 2 + c J m\ c 2 + (tt 12 - — 7T {12)3 ) 2 + 
V m V m 

H pr— + VdarwinKPu', ^12 H 7T(i2)3; ~^i2 H ^(12)3) = constant, 

4tt|p 12 | " m m 

2 

V ma tter(r) = ^ ^(r) 7T(i 2 )3 = constant, 
1=1 

Jmaiterli") = ^( r ) x ^( r ) = x ^12+ (v + l; P(i2) 3 ) x ^(12)3 = constant, 



(5.3) 



while from Eq.(4.2) we get 



2 y 

^matter = ~ ^ V m i ^ + *i ~ 



i=l 



IQ1Q2 

2 c 



ft 



I' {2 gfc 2A ± 52(«2,Pl 2 JJ 



77! * , - + 

' 2 2 , - 2 



+ V2 " 



m 2 c 2 + ft 2 



1 / / 2 2 , - 2 5 , 4/ 2 2 , - 2 9 f r — ^ 

- o ymfc 2 + «! — ^ + ym^c 2 + ft 2 — ^ /Ci 2 (ki,K2,Pi 2 ) - 

</ ^ la — ffo I la — 7T1I ' 



47T C 



Q1Q2 /" , 3 _ - / 



/ 



d a a (7r±si(^-7h,Ki) ■ n ±S2 (a - fj 2 ,FZ 2 ) + 



+ B si (a -rj^Kx) -B S2 (a-rj 2 , k 2 )) = 



— ^ ^matter Rc{ T ) ■ (5-4) 

C 



While the natural choice for the effective center of motion is center of energy ( c = R c = 
~ c fe^> Ec i-( 5 - 3 ) suggests the other possibility <f matter (r) =Vi2=V + \ P(i2)s wi th 77 given 
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in Eq.(4.11): it would imply the relation jf matter (r) = p 12 x 7r 12 + Cmatter ( r ) x ^(12)3 and , 

differently from the choice ( c = R c , (-matter is a canonical variable like ffu (however 7T(i2)3 is 
not the conjugate variable). 
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VI. THE DIPOLE APPROXIMATION OF ATOMIC PHYSICS 



Let us now look for the relativistic generalization of the electric dipole approximation of 
semi-relativistic atomic physics. To this end we consider a 2-particle system and we replace 
the particle 3-positions ffi and 3-momenta % = 1,2, with the naive center of mass and 
relative variables of Eqs. (2.3): fji 2 , k 12 , pi 2 , ^12- 

The standard electric dipole moment is directed along the relative variable pi 2 . Usually 
only neutral systems with opposite charges of the particles are considered in applications 
of the dipole approximation. As a consequence it is convenient to introduce the following 
notation for the Grassmann-valued charges 



Qi = Q + Q, Q2 = -Q + Q, Qi = Ql = 0, QiQ 2 = Q 2 Qi^0, 



Q = \(Qi-Q2), Q = \(Qi + Q2), Q 2 = -QiQ2, Q 2 = Q1Q2 if Q ^ 0. 

(6.1) 

The restriction to opposite charges is done by introducing the constraint Q ~ 0, which 
implies e = Q ~ Q\ ~ —Q 2 . This allows us to discard the terms in Qi Q 2 coming from Q 2 
from those coming from Q 2 : only these terms will produce effects of order e 2 in a neutral 
system with Grassmann regularization (it eliminates the e 2 terms coming from Q\ and Q 2 ,)- 

Strictly speaking it is only after the quantization of the Grassmann-valued charges Qi, 
sending each of them in a two- level system with charges (+e, 0) or (— e, 0), that we are really 
allowed to use e — Qi — — Q 2 . 

The natural definition of a Grassmann-valued electric dipole moment is 

d(r) = Qpi2(r) = < ^ 1 - ^ 2 pi 2 (r) ^ e=Ql ~_ Q2 ep 12 (r). (6.2) 



The alternative definition 



D(T) = ^Q i fH(T)=d(T) + 2Q[ff 13 (T)+ V i2(r)j ^ S ,o d(r), (6.3) 

i=l 

is equivalent to Eq.(6.2) for neutral systems. 

The electric dipole has not to be confused with the Dixon spin dipole, which is oriented 
along the spin S = pu x nu and determines the direction of the magnetic dipole moment. 

The length |pi2 (^~) | is of the size of the atom, i.e. a few Bohr radii. For a monochro- 
matic electromagnetic wave, i.e. for A±(t, 7712 (t)) ~ ae lk '^ 12<yT ^ with \k\ = we have 

ii x(T ,U))i 1 i Mr) ■ A) Mt > r/l2(r))l w pl2(r) ' ^ 27r """"fi^- This rati ^ is <K 1 in 

the /on^ wavelength approximation \ em >> size of atom (i.e. for radio-frequency, infrared, 
visible or ultraviolet radiation). 
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To get the dipole approximation [7] we use Eq.(2.3) and we make the following expansion 
(m 3 = mi) 

= A ± (r,ff 12 (r)) + (-Y+ 1 U 2 {r) ■ ■£-) MtMt)) + 

m V orivi' 

+ 0([p 12 -d ni2 ] 2 A ± ). (6.4) 
For the internal Poincare' energy generator of Eqs.(4.1) we have (// = m\ m 2 jm) 



£(int) = Mc 2 = 

= c m\ c 2 + k 12 (t) + tt 12 (t) j + ^/m 2 c 2 + k 12 (t) - tv 12 (t)^ ) 

_ ^ mi Qi ! m 2 Q 2 j 

ym?c 2 +(^/? 12 (r) + 7r 12 (r)) 2 ^/m 2 c 2 + (=* « 12 (r) - vr 12 (r)) * 

«12(t) 



A ± (r,7/i 2 (r)) 



m 2 c 2 + 2 12 (r) + vr 12 (r)) 2 ^/m 2 c 2 + * 12 (r) - tt 12 ( 
7Ti 2 (r) • A±(r,77i 2 (r)) - 



2 

r)" 
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Its semi-relativistic limit with the restriction e = Q\ ~ —Q 2 is 



p TOr .2 1 /&(t) ^12 (t) e 2 

t(mt) — >c->oo,e=Qi«Q 2 mC H o <~ 



2m 2fi 47t|pi2(t) 

e 7Ti 2 (r) 



c // 



^(r^i 2 (r)) 



e «i 2 (t) , v d 



( /7l2(r) '^;)^ (r ' r/l2(r)) - 



cm V dfji2 
e m 2 - mi tt 12 (t) /_ d 



(pi2(r)- 7 —)A ± (T,ff 12 (r)) + 



cm /i V <9f/i2- 

+ 0(c- 2 ) + 0([p 12 -d lJl2 ] 2 A ± ). 

(6.6) 

In this way we recover a semi-classical version of Eqs.(L3), (L4) and (14.34) of Ref. [7], 
without the e 2 terms corresponding to Q 2 = 0. 

If we could evaluate the internal 3-center of mass q + (see Subsection IIA) as a function 
of if 12 , K12 = k+, P12, 7Ti2, then the results of Section IV would allow us to write the internal 
Poincare' generators (4.1) in the following form (see Subsection IIA for the relative variables 



£( int ) = Mc 2 m c \J M 2 c 2 + k\ + {junction of q + ,K + ,p q , 7r 9 ) + 

+ v: m + o(p 12 2 ), 



V{int) = K + + V em ~ 0, 



Vl2 X «i2 + Pl2 X 7Ti 2 + J e[ 

— » 

— » — * . — * — * 

9+ X K.+ + Pg X 7T g + t/em) 



£(»nt) = /C[g + ,/t + ,pg,7Tg, ss 0. (6.7) 

^C(int) ~ determines q + — if + .. as it was done in Section IV for if. In this way we would 
get £{i n t) as a function only of the relative variables p q , n q in a form useful for the electric 
dipole approximation. 

While the external decoupled (canonical non-covariant) 4-center of mass P M , has an 

effective mass M = £( int )/c 2 and a spin S = Ji^nt), the particle subsystem (the atom) has 
the effective mass Ai of Eq.(6.7), a position q + , a 3-momentum k + and a spin given by 
the matter part of Ji^nt)- These quantities replace M c , ( c , V c and J c of the pole-dipole 
approximation to the multipolar expansion of the previous Section and give a canonical 
pole-dipole description of the atom. We have the following replacements 
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£(t) >-»• ?7i2(r) or <f+(r), 



"P c (r) «i 2 (r), 

^c(t) i—> Ai Q (relative variables) + Coulomb potential + interaction with the electro — 
magnetic field , 

Jc(r) h-> matter part of J^ int ). 

The results of Section IV allow us to eliminate the overall internal center of mass and 
introduce a dependence on the variables 7f(i2)3 = constant of motion and P( 12 ) 3 , describing 
the relative motion of the atom with respect to the collective variables of the electro-magnetic 
field configuration. 

Usually atoms are described not as the quantization of extended open subsystems with an 
effective 4-center of motion as in the previous Section, but as point-like systems (monopole 
approximation: M c i— > m, Cc(t) l— ► V12, 'Pc(t) | — > ^12) with some additional structure 
(higher multipoles) describing the energy levels and the interaction with an electro-magnetic 
field. After quantization this point-like description of positive-energy atoms leads to the 
effective Schroedinger equation of Ref. [14] used to describe the external propagation (its 
de Broglie wave) in atom interferometry. Conceptually this effective Schroedinger equation 
should be derived by studying the positive-energy sector of solution of some wave equation 
with a fixed mass and a spin (spin dipole) 18 , which couples to the magnetic field. As shown 
in Ref. [16], the description of charged positive-energy spinning particles in the rest-frame 
instant form can be made by using Grassmann variables for the spin. 

In this approximation with a fixed mass one looses all the internal structure of the atom, 
described by its energy levels. To remedy it one adds a finite-level structure to the point 
particle at the quantum level: this allows to consider more realistic approximations for the 
coupling to electric fields in the electric dipole approximation. The simplest model is the 
two-level atom approximation [7, 17]. 

As it will be shown in paper III, we can define a system in the rest-frame instant form, 
which after quantization leads to a two-level atom whose electric dipole interacts with the 
electric field after the transition to the electric dipole representation. There will be extra 
Grassmann degrees of freedom for the description of the two levels. 

Therefore, we must now find the relativistic generalization of the electric dipole represen- 
tation in the limit of the dipole approximation. 



For instance the Dirac equation is used for the coupling to external gravitational fields and the study of 
gravito-inertial effects [14]. 
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VII. THE RELATIVISTIC ELECTRIC DIPOLE REPRESENTATION 



In atomic physics the electric dipole approximation suggested the introduction of the 
electric dipole representation, where the interaction term A± ■ k { is replaced with the inter- 

— * — * 

action of the electric field with the electric dipole d of Eq.(6.2), d ■ rr±. In this way there 
is no explicit dependence on the transverse electro-magnetic potential: only the transverse 
electric and magnetic fields appear. 

This representation is a particular case of an equivalent formulation of electro-dynamics, 
as explained in Chapter IV of Ref[9]. Equivalent formulations are obtained 

a) by a change of the gauge of the electro-magnetic field (here we are using the radiation 
gauge, a special case of Coulomb gauge); 

b) by a unitary transformation e l s corresponding to a classical canonical transformation 
whose generating function S 

bi) is determined by a change of the particle Lagrangian by means of a total time deriva- 
tive with the electro- magnetic field considered as an external field with a given time 
dependence; the electric dipole representation is obtained with the G0ppert-Mayer uni- 
tary transformation (see p. 635 of the Appendix of Ref.[8] and pp. 266-275 of Chapter IV 
of Ref. [9]), which is useful when ^ Qi = and whose classical generating function is 
S = \ d{r) • A±(t, 7712 (t)) with d{r) of Eq.(6.2) (use is done of the dipole approximation) 

1!) 20. 

) 

62) is not determined by a change of Lagrangian but still with the electro-magnetic field 
considered as an external prescribed field 21 ; 

63) is connected to a change of Lagrangian with a dynamical electro-magnetic field but 
with the dipole approximation replaced by a description of the localized system of charges 
by polarization and magnetization densities 22 . 
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19 Strictly speaking the Appendix of Ref. [8] (see Eqs.(4) and subsequent one) makes the separation of the 
matter part from the radiation field in the Coulomb gauge and in the unitary transformation uses only 
the radiation electric field. Therefore our canonical transformation of I is the analogous separation in 
the radiation gauge (we get Darwin as an extra bonus!). The final Hamiltonian of Eq.(4.13) has the 
interaction reintroduced by the vanishing of the internal boosts (for the 1/c expansion see Eq.(4.5)). As 
a consequence a canonical transformation with generating function S = i d(r) ■ A± ra d( T , Vi2( T )) should 
now work like it happened in the Appendix of Ref. [8] . 

A generalization of the G0ppert-Mayer transformation is given in Ref. [9]. Its generating function is 
Sz = \ d(r) ■ £ dXA ± (r,Xff,(r) + (1 - A) rf 2 (r)) = \ d(r) ■ £ dX A ± (r, ft 2 (r) + p 12 ( T )) = 

\ d(r) ■ fl dX A ± (t, 7?i 2 (t)) - (^ - A) (pi 2 (r) • ^) A ± (t, tJ 12 {t)) + .. , with the integral taken along 
the straight-line joining the two charges. The first term is the G0ppert-Mayer generating function. 
See the Hennenberger unitary transformation, given at pp. 275-279 of Ref. [9], changing the positions 
instead of the momenta: the generating function is S = ^ Y^i ^ ^ ( T ) ' ^( T ' ^) w i tn Z{ T -> V) = 
— f T dr\ A ex ternai(Ti,ff). When the electromagnetic field is dynamical and not external, it is considered 
a quantized radiation field and the Hennenberg unitary transformation leads to the Pauli-Fierz-Kramers 
transformation. 

See the Power- Zienau-Woolley transformation of pp. 280-297 of Ref. [9] generalizing the G0ppert-Meyer 
one. 
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In Appendix L of Ref. [7], the semi-relativistic electric dipole representation is obtained by 
starting from the particle Hamiltonian of a neutral 2-particle system in the dipole approxi- 
mation (the electro-magnetic field is considered external), by evaluating the corresponding 
Lagrangian, by eliminating a suitable total time derivative ^ and by reverting to the Hamil- 
tonian. 



A. The Semi-Relativistic Electric Dipole Representation with Grassmann- Valued 
Charges 

In Subsection 1 of Appendix B the results of Ref. [7] on the electric dipole representation 
are revisited in the case of Grassmann-valued electric charges (Qf = 0, Qi Q 2 ^ 0) in 
the radiation gauge for the electro-magnetic field, which is treated as an external field, by 
starting from the Hamiltonian = = {F, H c}) 

H c a w)-*UrMTW = £ _ f a m . idT> m): m 

^ 2 mi ^ 2 mi ^ c m,, y ' /n ;; ' v ; 

i=i 1 i=i % i=i * 

As shown in Appendix B, the final Hamiltonian (B6) becomes the R0ntgen Hamiltonian 
(7.2) of the electric dipole representation, given in Eq. (14.37) of Ref. [7], for 2 Q = Qi + Q 2 ~ 
0. ( Q. 



+ 
+ 



12 



7' ) 71 



12 



2m 

— «12(t) 

mc 



2/i 



+ V(p 12 (r)) + 



pia(r) x S(r,?7 12 (r)) 



+ 



2 /x c 



7Ti 2 (r) • pi 2 (r) x £(r, 7712(7-)) 



m 



-P12W •vr ± (r,r7 12 (r)) - 

c 2c m 



mi 



3/i 
4m 2 



p 12 (r) x B(T,ff 12 (r)) 



(pn(r) • pi 2 (r) • 7r ± (r, ^(r)) - 



•95712 



(7.2) 



Note that the last term has not the coefficient ■£— 

in our calculation the terms Qj = are missing before imposing the condition Q 1 + Q 2 



of Eq.(15.37) or (L.14) of Ref. [7], because 

0. 



The resulting generating function S of Eq.(B3) is an extension to the next order of S Q , 
the generating function of the classical G0ppert- Mayer transformation. 



B. A Generating Function from a Relativistic Lagrangian with External Electro- 
Magnetic Field. 

Let us now try to define a relativistic electric dipole representation in the framework of the 
dipole approximation with the same method but starting from a particle Hamiltonian given 
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by i \£(int) — £emj in accord with Eq.(4.1). In Subsection 2 of Appendix B we determine 

a relativistic Lagrangian for the particles in the dipole approximation treating the electro- 
magnetic field as an external one, following the same method as in semi-relativistic atomic 
physics [7]. This identifies the generator of a canonical transformation, which will be used 
to find an electric dipole representation. 

If we use the dipole approximation for neutral systems, Q ~ 0, the emerging total time 
derivative identifies the following generating function S, given in Eq.(B14) and coinciding 
with S of Eq.(B3) at the lowest order, (Eqs.(2.3), (6.2) and (6.4) are used) 



+ 0([p 12 ■ d r1l2 ] 2 A ± ) = S\ Q *o + 0([p 12 • 4 2 ] 2 A ± ), 



■A 



»)] 



+ 



= ^i + o([p 12 -4 2 ] 2 i ± ), 



m 2 -mi d 



2 m dffi2 



p 12 {r) ■ A ± (r, 7712 (r)) 



Si 



2c 



d(r). A ± (r,fj 1 (r))+A ± (T,ff 2 (r)) + 0([p 12 ■ 4 2 ] 2 A ± ). 



(7.3) 



We also find that S] q~ an d ^Is^o differ by terms 0([p 12 ■ dfj 12 ] 2 A±) from a generating 
function 5i = j^d(r) ■ Y^=i A_l(t, ffi(r)), which could have been defined independently by 
using the dipole approximation. 

The inverse Legendre transformation from the resulting Lagrangian (B15) is the Hamil- 
tonian of Eq.(B20) with the Ai(r) given in Eq.(B16). If we add to it £ em , we get the 
Hamiltonian M e .d. r .c in the dipole approximation for the electric dipole representation re- 
placing V7 int -s = Mc of Eq.(4.1) of the standard representation (/%(r) are the new momenta 
of Eqs.(B17)) 
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M e . d . r C 



[Mr) 



£e.d.r. — 

C 

Ql-Q2 

2c 
Q1 + Q2 

2c L m 
m 2 — 2mi m 2 _ 

2^ Pl2 



~l v v \Jm\c 1 + Ki 2 (r)/ 



+ 



Q1Q2 



m 2 — mi 



2m 

Pi 2 (t) • 7r±(r,77i 2 (r)) 
9 



^Jm 2 c 2 + k1 2 [t)) 4:Tcc\p 12 (r)\ 

) • (pi2(r) • ^-7r±(r,77i 2 (r))) 



tt± (r, 7712 (r ) ) — — P12 (r) • 



/ a 
(r) • (pi 2 (r) • ^-7r±(r,7fi 2 (r)) 



• e ,o E (^^ + ^-v)- V m?c2 + ., 2(r) y , 4icfe(T) 



df) 12 

■ Air j \ 



If 

C 



^(r)-A(r) \ + QiQ 2 



Q1-Q2 
2c 



P12 (r) • tt± (r, 7712 (r) ) - ^ P12 (t) • (Vi 2 (t) • tt± (r, t7 12 (r ) ) ) + 



-)))] 



C 



Mr) ^s,o -^Pi2(r) x 5(r,r/ 12 (r)). (7.4) 

This Hamiltonian becomes the R0ntgen Hamiltonian ifi of Eq.(7.2) in the semi- 
relativistic limit. However it has not the form expected in the electric dipole representation. 

C. The Relativistic Electric Dipole Representations Induced by the Generating 
Functions S and S±. 

To try to find a relativistic electric dipole representation, let us study the classical 
point canonical transformation generated by T$ = e^' s ^ with the generating function S 
of Eqs.(B3)-(7.3) (coinciding with S in the dipole approximation) 23 by considering also 
the electro-magnetic field dynamical (Eq.(7.4) was obtained by considering it as an external 
field). In the canonical basis ifi 2 (r), ki 2 (t), /?l 2 (t), 7?i 2 (t), A±(t, a), 7c±(r,a), the variables 
Ai_, 7712, P12 remain unchanged. Instead for the momenta and the internal energy of Eq.(4.1) 
we get 



Since S — Qi S1+Q2 S 2 depends only on the coordinates, we get F = T s F = F+Qi {F, Si}+Q 2 {F, S 2 } + 
Qi Q2 [{{F,Si},S 2 } + {{F,S 2 },Si}] =F + Q 1 {F,S 1 } + Q 2 {F,S 2 }, since {{F,S a },S b } = 0. 
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- / \ t s -j , , — * / \ Q1-Q2 d 
k 12 {t) -> /c 12 (r) = «i 2 (r) - 



Pi 2 (t) • A ± (r,r/ 12 (r)) + 



m 2 - m\ ( _ , s (9 
+ — - Pi 2 (r 



2m 



2c 9 7712 

pi 2 (r) • A±(t,77i 2 (t)) 



7?i 2 (t) S 7T 12 (t) = 7ri 2 (r) - 



dff 12 

Q1-Q2 
2c 



A±(r,ff 12 (r)) + 



+ 



2m 



rrii 



^(r) = -^ 12 (r) + (-r 1 7r 12 (r), 
m 



_|_)i ±(T , ?12(r)) + _|_( &(r) .i ±(T , ?12(T)) )) 



ttKt, a) 5 ^(r, a) = n r ± (r, a) + p* 2 (r) P?(a) \s 3 (a - ff 12 (r)) + 



+ ^(p 12 (,)-4)^-^)) 



(7.5) 



If we apply T$ = e^' 5 ^ to £(i„ t ) of Eq.(4.1) we get its expression in the new canonical 
basis ?7i 2 (t), « 12 (r), Pi 2 (r), tt 12 (t), A±(r,a), i? ± (T,a) 



c °(int) 



2 

Af'c = V Jm 2 c 2 + zJHt) + - ^ . , + 
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i=i 
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kAt)-FAt) 
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<5i Q2 2 



2c 2 ^ v /m 2 c 2 + /5/ 2 (r) 

2c 
Q1Q2 



mf c 2 + k/ 2 (t) 
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[pi 2 (r) • n ± (r, ff 12 (r)) + ™" 1 (Vi 2 (t) • -^-j (pi 2 (r) • ttJt, 7712(1-))) - 



2c JM + 2c 



n ± ' 2 + B* 
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+ 



m dff 12 
m 2 — rri\ 



pi 2 (r) • A ± (T,ff 12 (T)) + 
d 



Pl2(T) 



2m V y (9 1712 
+ H l+1 [l ± (r,r/ 12 (r)) + 



m 2 — mi 



d 



p 12 (r) • A ± (r,77i 2 (r)) 



9 



+ 



I(t) = Um 3 _^ 2{T) p r 12 {T)p s l2 {T)Pl s {a) [<5 3 (a - rf 12 (r)) + 
m 2 — mi / _ , s (9 



+ (^) 2 (^)-4) 2 ^-- M) ]- 



(7.6) 



It can be checked that in the dipole approximation and in the semi-relativistic limit one 
recovers the R0ntgen Hamiltonian Hi (7.2) of the semi-relativistic electric dipole represen- 
tation: £' {int) = M' c 2 = m c 2 + H ± + £' em + 0(c- 2 ). 

However, at the order 0( ) a singular term 7(r) appears. It comes from the 7r 2 (r, a) 
term of the original electro-magnetic energy density. Notwithstanding the Grassmann reg- 
ularization, we get this diverging term as a consequence of the dipole approximation: both 
ffi(r) collapse in the collective variable r/i 2 (r) describing the motion of the electric dipole. 
Moreover the transverse vector potential is still present in Eq.(7.6) like in Eq.(7.4). 

The same problem arises if we consider the point canonical transformation T§ = e^' Sl ^ 

~ — * — * 

generated by the function Si — ^ d(r) ■ J2i ^J-( r ) Vi( T )) defined in Eq.(7.3) . Also in this 
case we get expressions of type of Eq.(7.6), which we omit, with a similar singular term. 

The use of S, instead of S or Si does not change the situation. Therefore it is not clear how 
to define a relativistic electric dipole representation, free from singularities and connected 
with the semi-relativistic one in the dipole approximation, when the electro-magnetic field 
is dynamical as it happens in the rest-frame instant form. 



D. A Relativistic Lagrangian with Dynamical Electro-Magnetic Field and the In- 
duced Relativistic Electric Dipole Representation 

However a relativistic representation in which the singular term is replaced by a contact 
interaction at order 0(c~ r ) can be defined by defining a Lagrangian for the isolated system 
"charged particles plus the electro-magnetic field" on the instantaneous Wigner 3-spaces. 
This is done in Subsection 3 of Appendix B. By modifying this Lagrangian with a total time 



In this case we have Si = (Qi — Q2) Si with Si = j^pi2{r) ■ J2i A±(T,rji(T)) depending only on the 
coordinates. Therefore we get F = T §i F = F + {Qi - Q 2 ) {F, Si} - Qi Q 2 {{F, Si}, Si}. 
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derivative, we can identify a new generating function S 2 = \ Y^h=i QiVi( T ) ' Ai_( T > Vi{ T ))i 
whose dipole approximation is given in Eq.(B28), leading to this new representation. 

Again the variables A±(r, a), rf 12 (r), pi 2 (r) remain unchanged under the point canonical 
transformation generated by T$ 2 = e^' s ' 2 ^ with generating function S 2 - Instead for the 
momenta we get 25 



nl(r,a) ^ ^(r, B) = tt^t, a) + £ Q, P?(a) S 3 (a - #(r)) i£(r), 

i=i 

«r) " f A 1( r, *<r» 5» * (r) + ^ *<r) • 



rf 3 o" 7r]_ (r, <r) — > / <i 3 (T 7r^_ (r, a) — 



2 £ g.r/^r) -^(r,^(r)) + 2g 1 g 2 r ? [(r)^(r)Pr(7/ 1 (r))5 3 (7/ 1 (r) -^ 2 (r)). 
i=i 

(7.7) 

The internal energy of Eq. (4.1) is replaced by the following expression 
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We 



use F sr = -e s ™ B u , R\ V f ^ = «f r?| (^ sr + ^f ) = x B+k? = -ft x ^-B+Sf r/f 
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(7.8) 



The Lorentz boost K,( int ) of Eqs.(4.1) is replaced by the following expression 
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(7.9) 



The new form of P(j„t) and J(mt) of Eq.(4.1) after the canonical transformation is 
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T) X Ki T + 



i=l 



i/«<Wx( 



jfixB (r,(j 



(7.10) 



In Eqs.(7.8)-(7.10) the transverse potential Aj_(t, a) has to be re-expressed in terms of 
the magnetic field B(t, a). 

In Abelian electro-magnetism, where the 2-form F rs da r da s is closed, the vector potential 
Ap(r,a) of the radial (or Poincare') gauge, satisfying a ■ Ap(r,a) = 0, leads to the exact 
1-form Ap{r, a) • da, so that the Poincare' lemma implies 



Ap^fftir)) = - / dXXf)i(T) x B(T,\ffi(T)). 



(7.11) 



The transverse potential Aj_(t, a) of the radiation gauge can be connected to Ap(r, a) by 
means of a gauge transformation such that A±(t, a) = A p (t, a)+d A(r, a) with the following 
consequence of transversality 



AA(t,<t) = -d-A P (T,a) 



d x B(t, A a) 



= - j d\ A a ■ 

J o 

whose solution is (we disregard solutions of the homogeneous equation) 

A(t, a) = J dVi 



lo d\\(Ti- 


— * — * 

<9i x B(t, A(?i) 
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(7 — <7i| 



Therefore we get 



(7.12) 



(7.13) 



A ± (r,a) = - / dXXa x B(r,Xa) + 



+ d / dVi 



Jo 1 dX X ff i 



d l xB 



(r,A^)] 



47T <7 — (jd 



(7.14) 



The term depending on A_i_(t, a) in the energy (7.8) and in the boost (7.9) becomes (we 
use the spin notation Si = ffi x Ki) 
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(7.15) 



With this equation we can also evaluate Ki(r) ■ (rji{r) ■ A±(r,ffi(r)) = f)i(r) ■ (^i(r) ■ 
4)l ± (r,^(r)) + S ? l (r).fi(r,^(r)). 

As a consequence we get a generalized electric dipole representation: besides interactions 
of the particles with the transverse electric field at their position and interactions with a 
delocalized magnetic field ( one of them is with the particle angular momentum), there 
is a contact 0(c _1 ) term replacing the divergent term arising from the previous canonical 
transformations. 

If we make the dipole approximation and the semi-relativistic limit of S^ nt ) we get 
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(7.17) 



With respect to the R0ntgen Hamiltonian (7.2) there are numerical factors \ of difference 
and extra terms containing the transverse vector potential and effects due to the motion of 
the collective variable f/i 2 (r). 

Eq.(7.8) seems to be the simplest relativistic extension of the semi-relativistic electric 
dipole representation in the framework of the rest-frame instant form. 
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VIII. CONCLUSIONS 



In this second paper we concluded the classical treatment of the isolated system of 
charged positive-energy scalar particles with Grassmann-valued electric charges with mu- 
tual Coulomb interaction plus the transverse electromagnetic field in the radiation gauge in 
the rest-frame instant form of dynamics. In the first paper we showed how every isolated 
system (particles, fields, strings, fluids) is described in the inertial frame centered on its 
Fokker-Planck external center of inertia: the isolated systems is replaced by a decoupled 
canonical non-covariant Newton- Wigner center of mass (described by the frozen Jacobi data 

— * 

z, h) carrying its internal mass and spin in such a way that there is an associated external 
realization of the Poincare' group. The dynamics inside the isolated system is defined by 
construction in the instantaneous Wigner 3-spaces and is described only by means of relative 
canonical variables. This happens because the internal realization of the Poincare' group 
built from the internal energy momentum tensor of the system is unfaithful: there are the 
rest frame conditions V{i nt ) ~ and its gauge fixing JC( int ) ~ (vanishing of the internal 
Lorentz boosts), whose role is the elimination of the internal 3-center of mass and of its 
three momentum. 

In this paper we have studied how to find collective variables in the instantaneous Wigner 
3-spaces such that it is possible to solve the equations K,^ n t) ~ explicitly. While for particle 
systems there are no problems, the construction of collective variables for the configurations 
of massive and massless fields is nontrivial and only holds for a subset of the field config- 
urations with finite Poincare' generators. We adapted to the rest-frame instant form the 
methods developed by Longhi and collaborators in Refs.[5, 6] and used them to find a so- 
lution of the equations JC( int ) ~ for the transverse radiation field. This made possible to 
solve this equation also for the isolated system of charged particles with mutual Coulomb 
interaction plus the transverse electro-magnetic field after the canonical transformation of 
paper I. 

The final result in an arbitrary inertial frame is that inside each instantaneous Wigner 
3-space, with origin in the external covariant Fokker-Pryce center of inertia, the internal 
3-center of mass becomes a function of the relative canonical variables, so that the dynamics 
is described only in terms of them. 

If we consider an isolated 2-particle system (but the same conclusions hold for an arbitrary 
isolated system) we have the following situation on the instantaneous Wigner 3-space: 

A) the particles are described by the two Wigner 3- vectors fji(r) = fji 2 (t) + 
R m £pi 2 (r) « f(pi2(r),n 12 (r)) + ^p 12 (r) with p 12 (r) = ft(r) - ff 2 (r); 

B) as shown by Eq.(I.2.19) the external non-covariant canonical center of mass, a? M (r) = 
Y^(t) + e£ (h) a r , is identified by the 3-vector a = sxh ^y e cou y ^hen try to define 

two new particle 3- vectors ff^r) = a+ (— ) i+1 ^ Pui^) and describe the dynamics with them 

and their conjugate momenta rather than with z, h, pi 2 , ~K\ 2 : this would be equivalent to 
the non-relativistic treatment with x< n \i and pi n u. However the obstruction to do this is 
the non-covariance of the external canonical center of mass, which implies that a is not a 
Wigner spin-1 3-vector. The presence of interactions in the Lorentz boosts and the associated 
No-Interaction theorem, which imply the existence of three relativistic collective variables 
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replacing the unique non-relativistic center of mass and spanning the non-covariance M0ller 
world-tube around the external Fokker-Pryce center of inertia, lead to a description in which 
the decoupled canonical center of mass breaks manifest Lorentz covariance and is a global 
quantity which cannot be locally determined. 

Then we clarified how to get the multipolar expansion of the energy-momentum of the 
particle subsystem and its non-canonical pole-dipole approximation, not to be confused with 
the canonical pole-dipole structure carried by the external Newton- Wigner center of mass. 
Thinking to strongly bound clusters of particles as classical atoms, it is found how to define 
their effective 4-center of motion. 

Finally we have studied the relativistic electric dipole approximation. We have shown that 
it is non trivial to find a relativistic generalization of the electric dipole representation used 
in atomic physics. When the electro-magnetic field is dynamical, as requested by the rest- 
frame instant form, the point canonical transformations suggested by the semi-relativistic 
approach tend to generate singular terms at the order 0(c~ v ). They are connected to the 
use of the dipole approximation. From a study of the Lagrangian on the instantaneous 
Wigner 3-spaces it is possible to define a relativistic representation in which the singular 
terms are replaced by contact interactions. In this new relativistic representation, whose 
semi-relativistic limit is different from the standard electric dipole representation, there are 
interactions only with the transverse electric and magnetic fields, like in the semi-relativistic 
case. 

In the third paper we will quantize the rest-frame instant form of our formulation of 
relativistic particle mechanics. This will allow us to consider relativistic bound states of 
clusters of particles as atoms. The coupling to a quantized transverse electro-magnetic fields 
will allow us to define relativistic atomic physics and relativistic generalizations of the two- 
level atom. Future applications will be relevant to study relativistic entanglement and atom 
interferometry. 
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APPENDIX A: DEFINITION OF THE COLLECTIVE VARIABLES FOR THE 
MASSIVE KLEIN-GORDON FIELD. 



The first formulation of collective and relative variables for a field was done in Ref.[5] 
for the free massive Klein-Gordon (KG) field. It was then adapted to the rest-frame instant 
form in Ref.[13]. Here we give the main results to introduce the methodology to be used to 
find such variables for the transverse radiation field starting from Ref.[6], where the results 
of Ref . [5] were extended to a scalar massless field. 



1. The Massive Scalar Field 

In Ref. [13] the scalar massive field (j)(x) was reformulated in the rest-frame instant form 
as the field <f>(r, a) = <P{zw{t, a)) on the Wigner hyper-planes z^(t, a) = Y^{t) + e£f(/i) a r . 
Its conjugate momentum is 7r(r, a), {0(r, a), 7r(r, ffi)} = 5 3 (a — ffi). A first canonical trans- 
formation replaces these canonical variables with their Fourier coefficients a(r,k), a*(r,k), 
whose non-vanishing Poisson brackets are {a(r, k), a*(r, k{)} = —i{2n) 3 2u)(k)8 3 {k — ki), 

uj{k) = \J m 2 c 2 + k 2 , k = 26 . 

Then a second canonical transformation replaces the Fourier coefficients with the 

— » — * — * — * / j*\ 

modulus-phase canonical variables I(r,k) = a*(r,k) a(r,k), (f(r,k) = jr In i , whose 

non-vanishing Poisson brackets are {/(r, k), </?(r, k-i)} = (2ir) 3 2 uo(k) 5 3 (k — k{). The final 
expression of the original fields is {dk = ^^^{k) ) 

0(r,a) = J dk \J / /(r, k) (^^,k)-i[^k)r-k-8\ + e -i<p(r,k)+i[u{k)T-k-8]^ 

7t(t, (j) = -i J dku(k) \Jl{T, k) ^ e iv(r,k)-i[u>(k)T-k-ff] _ e -i V (r,k)+iMk)r-k-ff]y ^ 

The internal Poincare' generators are (we must have a\(r, k), da\(r,k) e L 2 (R 3 ,d 3 k) 
for the existence of the following ten integrals ) 

VI = M 4> c = J dku(k) I{r,k), = J dkkl(r,k)^0, 

J<t> = = / dk I(r,k)k x -^<f(r,k), 
J dk 

/C = -V^t- j dkcu(k) I(r,k) -^<p(r,k) « 0. (A2) 
J d k 



The quantity Ua <J A = uj(k) r — k ■ a is built with o A = (r; (?) and k A = (fc T = uj(k); k). It is a Lorentz 
scalar because on the instantaneous Wigner 3-spaces the 3- vectors a and k are Wigner spin-1 3- vectors 
and k T = Lo(k) is a Lorentz scalar. For the free KG field the Fourier coefficients and modulus-phase 
variables are r- independent: in this Appendix we leave a formal T-dependence, because the formalism can 
be used also when interactions are present as in Ref. [13]. 
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In particular we have 

Mr,k),V A } = -k A , (A3) 

where V A = (7^; V*) and k A = (u(k); k). 

The collective canonical variables are assumed to have the form 

X; = J d~kcu(k)F T (k)p(T,k), Vl = M*c, 

X = J dkkF(k)cp(r,k), V^O, (A4) 
where the two Lorentz-scalar functions F T (k), F(k) have to be determined by the canon- 

icai conditions {x;,p;} = -l, {x;,v$ = 5 rs , {x;,x;} = {x;, x;} = {x;,v r ^} = 

{X^V^} = 0. The solution of these conditions is [13] 

mc A; 2 Vm 2 c 2 + k 2 mc k A 

J dku 2 {k)F T {k) = 1, y dku(k)k r F T (k) = 0, 

y dfc F fc s = -<T S , y dA;cj(A;) fc r F(ifc) = 0. (A5) 



This result implements the original construction of Ref. [5] adapted to the rest-frame in- 
stant form. The basic idea is to consider a real Lorentz-scalar function J-{k, Vf) normalized 
in the following way 

Jd~kk A F(k,Vf), k A =(u(k);k)=V A . (A6) 
Then, by using Eq.(A3) and (A6), one can check that the collective variable 

X^ = fdk dF ^K (r,k), (A7) 

satisfies {X},X»}=0, {X},Vf} = V B . 

Eqs.(A5) imply that the function F{k, Vf) has the form 

F(k, V*) = u(k) F T (Jz) V; - F(k) k ■ V*. (A8) 
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Since the functions F T (k) and F(k) are singular at fc — 0, the collective variables exist 

— * 

for field configurations whose phase variables satisfy </?(r, k) ^k^o k v with rj > 0. 

Let us remark that for field configurations 0(r, <?) such that the Fourier transform 

^ — * — * 

0(r, k) has compact support in a sphere centered at /c = of volume V, we get = 

The canonical relative variables, having vanishing Poisson brackets with the collective 
ones, are [5, 13] 



H(t, k) = dqg(k,q) I(r,q) - u(q) F T (q) / d~q x u( qi ) I(r, q x ) + 



+ 



F(q)q- J dqx q ± I(t, (fi)] , 



K(r,fc) = Vrip(r,k), 



{H(r, fc), K(r, fci)} = (2tt) 3 2 w(fc) <5 3 (fc - h), 
{H(r, fc), H(r, &)} = {K(r, fc), K(r, &)} = 0. 



(A9) 



Here (see Ref. [5]) G(q,k) is the Green function of the operator = 3 — m 2 c 2 A LB = 
3 - m 2 c 2 dl -2q- fy- (q-dj) 2 21 : V 9 Q(q, fc) = (2tt) 3 2 w(fc) 5 3 (fc - q). 

Another needed distribution is A(fc, q) = (2tt) 3 2 w(fc) 5 3 (k-q)-F T (k) u(k) u(q)+F(k) k-q 
enjoying the semi-group property J dq A(fc, <f) A(g, fci) = A(fc, fci). 

The final expression of the original fields in terms of the collective and relative canonical 
variables is 

<p(r,a) = 2 J dkA(T,k)cos(k-a + B(T,kj^, 
7t(t,ct) = -2 J dku(k)A(T,k)sin(k-a + B(T,kj^, 

A(r,k) = y]l{r, k) = sjF^k) co(k) V% - F{k)k-V^ + %H(r, fc), 
B(t, k) = (p(r, k) — u(k) r = 

= -h-Xt- u(k) (r -X;) + J dh dk 2 K(r, h) g(k u fc 2 ) A(fc 1? fc). 

(A10) 



27 Alb is the Laplacc-Beltrami operator of the mass-shell sub-manifold of the KG field. The operator 
T>^ is scalar, formally self-adjoint with respect to the Lorentz invariant measure for the massive case and 
has k A as null eigenvectors, V^k A = 0. 
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The classical field configurations which can be described in this way must belong to a 
function space such that: 

a) the internal Poincare' generators (A2) are finite (this is required in every approach to 
be able to build the Poincare' representation); 

b) the collective and relative variables are well defined: as a consequence of the existence 
of the Poincare' generators and of the study of the zero modes of the operator done 
in Ref.[5], we must have: b\) |7(r, k)\ — > k~ 3 ~ a , a > 0, \ip(T,k)\ — > k for k — > oo; b 2 ) 
\I(t, k)\ -> Ar 3+e , e > 0, \ip(r, k)\ -> A;'', r] > 0, for fc -> 0. 

The field configurations satisfying these conditions may be described in terms of the 
following multipolar structure in the rest-frame instant form: 

A) An effective particle X^ = (X^; X^) whose conjugate momentum = (PJ = 
Mfi c; ~ 0) describes the mass and momentum components of the Dixon monopole [13] 
in the multi-polar analysis of the energy-momentum tensor of the KG field. While X^ is 
a 3-center of phase for the field configuration, X^ is an internal time variable conjugate to 

= M^c. For particle systems an analogue of X^ does not exist, because the invariant 
mass is a function of the particle canonical variables and not an independent function like 
M^. As shown in Ref. [13], the rest-frame Hamilton equations imply that = M^c is a 

constant of motion and that X^ satisfies the equation — 1- 

B) A set of canonical relative variables H(r, k), K(r, k), which are constant of motion 

for a free field due to the complete integrability (Liouville theorem) of the free KG field. In 

terms of them we can build canonical multipoles with respect to the Fokker-Pryce center of 

inertia, origin of the instantaneous Wigner 3-spaces, and then all Dixon multipoles [13]. 

— * — » — * 

C) The canonical variables X^, V<j, are not independent: we have ~ from the rest- 
frame conditions and X^ — J dkuj(k) H(r, k) ^ K(r, k) from the elimination of the 

— * 

internal 3-center of mass implied by JC<f> ~ 0. 

D) The variables = M^c and X^ have the following interpretation. If we consider 
a constant energy surface Se, i.e. M^c 2 ~ E, in the KG phase space, then £ E is not 
a symplectic submanifold due to the presence of the extra variable X^. However, if we 

add the r-dependent gauge fixing X^ — a t rs (AJ= — 2 r from Hamilton equations) 
to the constraint c 2 — E m and we go to Dirac brackets, then Se becomes a phase 
space spanned by the canonical variables H(r, k), K(r, k) for every value of E. The Dirac 
Hamiltonian becomes a numerical constant Hp = E/c and, as shown from Eqs.(A6), the 
fields become r-independent (they only depend upon constants of motion for the free KG 
theory) . 

2. The Massless Scalar Field 

In Ref. [6] it was shown that the construction based on Eqs.(A6) and (A7) also holds for 
the massless scalar field. Therefore we have to find the limit of the function !F{k, Vf) for 

— * 

m — > 0. Now we have u(k) = k = \k\. 
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By using the definition of the delta function in the form S(x) = lim e ^ (^7= e x2 / e j 

with e = m 2 c 2 /Atc — > 0, we get that the functions F T (k) and F(k) of Eqs.(A5) become the 
following distribution in the massless limit 

8 7T 2 24 IT 2 

F T {k) = -gj- 5(k), F(k) = -— 5(k). (All) 

It can be checked that the conditions contained in Eqs.(A5) are still satisfied since 5(h) 
is an even function of k and we have 

/poo poo 
dkuj 2 (k)F T (k) =2 dkS(k)= / dk6(k) = l, 
JO J -00 

/poo poo 
dk k r k s F{k) = -2 5 rs / dk 5(k) = -5 rs / dk 5{k) = -5 rs . (A12) 
J0 J -00 

The collective variables (A7) are now well defined for the field configurations whose phase 
variables behave as </?(r, k) -^k~>o k and have the form 

1 r d^k 

S r d z k -> 

— J -^-k6(k)tp(r,k). (A13) 



Xa, — 



As shown in Ref. [6] the relative canonical variables, having zero Poisson brackets with 
the collective ones, are again given by Eqs.(A9), but now uj{k) = k and Q(k, q) is the Green 
function of the operator V$ = 3 — 2 q ■ — (q ■ d ( 



\2 28 



Also Eqs.(AlO) for the reconstruction of the massless KG fields hold. 

Now, the field configurations having well defined collective and relative canonical variables 
are restricted by the requirement of the existence of the Poincare' generators. This requires 
[6]: i) \I(r,k)\ -> AT 3 - 5 , 5 > 0, and \ip(r,k)\ -> k for k -> 00; ii) \I(r,k)\ -> k~ 2+e , 

— * 

e > 0, and \(p(r, k)\ — > k a , a > —5 for A; — > 0. However an extra infinite set of integral 

— * — * 

restrictions comes from the study of the zero modes Wi m (k) of the operator done in 
Ref. [6]. If we define the quantities Pi m = J dkwi m (k) I(r, k), then the restrictions are 

Pi m = j dkwi m (k)J-'(k,V^) for I > 2 (the Pi m with / > 2 are called super-translations, 
being connected with the generators of the algebra of BMS group studied in this framework 
in Ref. [15]). 



28 Like in the massive case this operator is scalar, formally self-adjoint with respect the Lorentz invariant 
measure for the massless case and has k A = (\k\; k) as null eigenvectors. 
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APPENDIX B: RELATIVISTIC LAGRANGIANS FOR THE ELECTRIC 
DIPOLE APPROXIMATION 



1. The Standard Electric Dipole Representation with Grassmann- valued Electric 
Charges 

Let us revisit the standard electric dipole representation starting from the Hamiltonian 
(7.1) with Grassmann- valued charges, we go to the collective and relative variables (2.3) and 



we make the dipole approximation (6.2) (i.e. we neglect terms of order Oy[p\2 • 
we get from the Hamilton equations (p = mi m 2 /m, a = r = ct, J dr L(t) = J dt L(t)) 



d \2 



A 



k 12 (t) = mcff 12 (r) + — (pi 2 (t) 

c V 



d 



drfi 2 



Q ™2 -m 1 / d \ - ^ 

c m V oi] 12 / 

vfi 2 (r) = /icp 12 (r) + ^ \A ± (r,ff 12 (r)) + — (pi 2 (r) • 

C L III \ U l\\2 



A 



+ 



+ « 

c 



m 2 — mi 



+ 77T.2 



m 



77i 2 (r)) + 

m z 

where the notation of Eq.(6.1) has been used. 



(Mt)~)a±(t,Vi2(t))], (Bl) 



dff 12 



Then by means of the inverse Legendre transformation we get the Lagrangian 



mc -J tic -J2 

-YVi2( t ) + y Pl2^) + 



+ i- 

c 



171/ \ or\\2 

2 2; n t- ( _ ( ** (Q Q m 2 -m 1 \ u 

c \ c c m 

Q m 2 — mi Q ml + m 2 \ u , s /_ , s 9 



m 



+ 2 % ) Pi 2 ( r ) • (Pi2(t) 



)p\ 2 (r)-l ± (r,^ 12 (r)) + 
AL(r,?7i2(r)). 



(B2) 



Let us write L = L\ + ^ with the following function 5 



5 = 
+ 



m 2 Qi - m 1 Q 2 

2mc 
m\ Qi+m\Q 



p 12 (r) ■ A ± (r,ff 12 (T)) + 



2m 2 



2 Pn(r) ■ ( Pn(r) ■ t^t- ) A ± (r, ff 12 (r ) ) 



•9^12 



Q 

c 



-* / \ r f - f . m 2 -m 1 _ 

pi 2 (r) • j4l(t, 7712 (r)) + — p i2 

2m 



( r ) • (V12O 



9 



+ — 



Q \ m 2 ~ m 1 _ -> _ mj + m 2 2 _ 

Pi2[T) ■ A ± (r,n 12 (r)) + 



m 



drfu 

2m- P^-{P^> dffl2 



AAr,ff 12 (r)) 
d 



+ 



>q*o — Mr) ■ \a±(t, ff 12 (r)) + m2 mi (pi 2 (r) ■ -J^-) A ± (r, ff 12 (r)) 
c L 2m \ omo/ 



) ■ (pi 2 (r) • A ± (r,77i 2 (r)) 

(B3) 
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where we used the notation of Eq.(6.1). For neutral systems, Q « 0, S — S Q + o(^(pi 2 



q4^) A±J is the extension to the next order of the classical counterpart S a of the generator 
of the G0ppert-Mayer unitary transformation. 



If we use dA ±{T ^(r)) = -jf^rMr)) + (fj 12 {r) 
lowing form 



A±(t, 7712 (t)), Li takes the fol- 



mc -j. . . fic 3. . . 2 Q ^ r _ 

£1 = -y^Or) + yPi 2 ( r ) + — ^12^) • Al(t,»7i2(t)) + 

'Q , Qm 2 -mi 



+ - + 

V c c m 



Pi 2 (r) -ttiXt, 7712(1-)) -rTi 2 ( r ) ' (/foCr) x B(t, rj 12 (r)) ) 



+ 



'Q m 2 -m l Qm\ + m|\ r / 9 \ ^ 

,7 -^r - + 7 -^H h (r) • l Pi2(r) • ^) ^ (r ' ^ i2(r)) - 



- Pi 2 ( r ) • (Pi2(r) X B(T,ff 12 (T)) 



+ 



(B4) 



We have used 5(r, ry 12 (r)) — ^4^- x A ± (r,77 12 (r)), p 12 (r) x B(t^ 12 (t)) = [pi 2 {t) 
gf-) -4±(T,?7i2(r)) - (Vi2(r) • Aj_ (r, 77i 2 (t)) ) . The new momenta are 

k 12 (t) =mcrf 12 (T) + —A ± (T,ff 12 (r)) - ( — + — m2 - mi ) p 12 (r) x B(r,f] 12 (r)), 



c c 

2 1 ™2 



-, . , Q m 2 - mi Q m{ + m 
n 12 (r) = »c Pl2 {r) - ( - + - 



pi 2 (r) x B(r,ff 12 (T)). 



(B5) 



Then by Legendre transformation we get the Hamiltonian of the electric dipole represen- 
tation in the dipole approximation 



#1 



K 12 [T) 7T 10 (T 



12 V 



2m c 



2/i c 



— + P12 r) • tt ± (r, 7712 (r) ) 

c c m 



+ 7 2m 2 J P12 ( r ) • (P12 (t) • ) t± (t, »7i2 (t) ) + 



c 2m 

1 /Q Q m 2 - mi \ 



+ — - + 

mc V c c 



m 



1 (Qm 2 — m\ Q m? + m 2 
' _| 



/i c V c 2m 
2Q 



•95712- 

) «i 2 (r) • (pi 2 (r) x B(t, 7712(1-))) + 
r 2//;2 ) ^i2(r) • (pi 2 (r) x B(r,77i 2 (r))) 



mc z 



K12W ■ i±(T,if 12 (T)) + (g+Q m2 m "" ) I±(t,t? 12 (t)) ■ (p 12 (r) x B(T,^ 2 (r))) 



For 2 Q = Qi + Q2 ~ 0, e = Qi, we get the R0ntgen Hamiltonian (7.2) of the electric 
dipole representation, given in Eq. (14.37) of Ref.[7]. 



61 



2. A Lagrangian for the Internal Motion of the Particles with an External Electro- 
Magnetic Field 



Let us now follow the method of Appendix L of Ref.[7] of determining the Lagrangian 
from the Hamiltonian by considering the electro-magnetic field as an external field. 

We begin with the internal energy \£(i n t) = 'Pjint) = Mc of Eqs.(4.1) and we do the 
inverse Legendre transformation on the particle variables to find an effective Lagrangian 
for the particles. After the dipole approximation and by using e = Q\ ~ —Q 2 , a term in 
the resulting Lagrangian will be a total time derivative identifying the generating function 
S. Then we will do the Legendre transformation of the Lagrangian without the total time 
derivative. The semi-relativistic limit of the resulting Hamiltonian turns out to be the 
R0ntgen Hamiltonian (7.2) of Ref.[7]. 

By ignoring the field energy £ em in Eq.(4.1), the inverse Legendre transformation produces 
the Lagrangian 

T ( \ - dff^r) dff 2 (r) 

= -TOicWl- I ""V - 1 -m 2 c\\\- I " '7 ' I , J?\^ 2 1 




4ttc|57i(t) -7f 2 (^) 

+ — A ± (t, 771 (r)) • — + — A ± {t, t] 2 {t)) ■ — . (B7) 

Switching to the center of mass and relative variables of Eqs.(2.3) (here we put fj = fj 12 
and p — P12) we obtain 



dfftr) m 2 dp(r) \ 2 f dff(r) mi dp(r) x 



' \ (I; 111 (li J 2 V V (h hi (It 

Qi Qi 



An c \p{r) 

Qi ( t i \\ , m 2 d x \ f dff(r) m 2 dp(r)\ 

+ 9l { A ± (t, m - ^ P(t) ■ A ± (r, fUr))) ■ (^1 - ^ *M ) , (B8) 
c m or] \ clt m ar J 

in which we have expanded the vector potentials about fj by using the dipole approximation. 
By using the notation Q — \ (Qi — Q 2 ), Q — \ (Qi + Q 2 ), of Eq.(6.1) after some algebra we 
find 29 



29 Note that wc divided the ■ (J){t) ■ A±(t, ?f(r))^ term into separate and equal parts following the 
treatment of Ref.[7]: this will produce the correct \ factors. 
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L(t) = — mi c y 1 

Q1Q2 
4nc \p(t) 



dfftr) m 2 dp(r) 
dr 



m dr 



+ 



m 2 c \ I i 
Am 



drftr) mi dp{r) 



dr 



m dr 



+ 



dfj{r) 



+ 



dr 
Am dp(r 



— (p- A ± (T,f)(r))^ - p(r) • —A ± (r,f)(r)) 



2m dr 

Q 



+ — • n ± (r, fj(r)) + p{r) ■ — (p ■ 7f ± (r, ff(r)fj - 

d_ 

dff 

d r , ^ „ d 



+ 



p(r) ■ —A ± (T,ff(r)) - g^[p(r) • A ± (r,^(r)) 
Am 



- — [fftr) ■ tt±(t, fj(r)) + p(r) • tt ± (t, ff(r)) + 

,2 



+ 



+ 



m — 2pm _ 



2m 2 



m 2 — 2pm dp(r) 



2m 2 



dr 



p(r) ■ ^A ± (r,ff(r)) - — (p(r) • A ± (r,ff(r)fj 



+ 



Am dp(r) 



m rfr 



P(r) • ^.4 ± (r,,7(r)) - A (p( T ) . .| ±(r , , 7 ( T )))' - 



A ± (r,77(r)) + 



dff 
dS_ 
dr~ 



(B9) 



in which we have ignored the second order spatial gradient terms of the vector potential 
(dipole approximation) and we have put Am = m 2 — mi, p = mi ™ 2 ■ 

The total time derivative terms, not contributing to the equations of motion, identify the 
generating function S: 



S = 



Q 



/\ 777 / Q _> \ " 

p(r) ■ A ± (t, ff(r)) + — p(r) • (p(r) • — A ± (r, ff(r))j 



c 



2m 
Am 



dff 



+ 



ff(r) ■ A ± (t, ff(r)) + — - p(r) • A ± (r, jftr)) + 



+ 



m 2 — 2pm _ 



2m 2 



— * 

Since we are in the radiation gauge we have E±(T,ff(r)) 



dA ± (T,ff(T)) 

dr 



(BIO) 
. Then by using 
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p x B 
rfx B 



dr 



d 



P- A ± (r,ff) 



P "dfj A±(T '^' 



7/ 



dfj 

dfj / d 
dr \dfj 

d d - 

A±(r,fj) - —fj- A ± (r,fj) -fj- — A ± (r,ff), 



A ± (r,fj))-A ± (r,fj)), 



<9r/ 



(Bll) 



we obtain 



,'dr/(r) m 2 cip(r 
Li(r) = -m 1 c\ 1 - I — ; 1 



Q1Q2 



Anx c |p(r 
dfftr) 



+ 



dr m dr 
Q s ^ , „ A/// 



m 2 c\ 1 



dff(r) vrt\ dp(r) 



— (p(r)-n ± (T,rf(T)) + T^P( T ) ■ (V( r ) ' ■g^^±(r,v(r)) S j 



dr m dr 
d_ 
dfj 



dr 
Q /Am 



p(r) x B(r, r7(r)) _ ^ Mil . fl T ) x b(t, fj{r)) ) + 



2m dr 
m 2 — 2pm 



,2 



• (p( 



<9 



(9 17 



*)))- 



m — 2pm dp(r) 

2m 2 dr 
Am dfjij) 



p(r) x B(r,fj(r)) - 



p(r)xB(r,fj(r)) + 



m dr . ^ r 

In the case e = Q Qi ~ — Q2, <2 ~ 0, this becomes 



(r,^(r))). 



(B12) 



^i(r) 



—mi C \ 1 

,2 



d?y(r) «i2 dp(r) 
dr m dr 



-m 2 c\ 1 



dfjir) m\ dp{r) 



+ 



r—— + -p r • ir ± (r,?7 r + - — 

47rc |p(r)| c c 2m 



dr m dr 
P( r ) ' (V( r ) ' ^^±(r,ff(r))^ - 



+ 



e dfjir) . e Am dp(r) _. . ->, _. . . 

p(r) x B(r,ff(r)) - - — - • p(r) x B(r,rf(r)). (B13) 



c dr c 2m dr 

and the generating function S reduces to 



5 = -p(r).[^(r^(r))+— 



9 



,<>)• — )i ± (r,77(r))) 



(B14) 



In order to obtain the Hamiltonian from Eq.(B12) we revert from the composite variables 
back to the constituent variables, which mean that we need to rewrite L in terms of the 
individual velocities using fj — fj 12 — miVl+ni2 m ; p — p 12 — f)± — ff 2 . 
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We obtain 



I,(r) = ^l-(^)^^1-(^V + 

Q1Q2 , Q ^ , . , ^ , N , Q Am _ . 9 _ 



+ — p(r)-7r ± (T,f](r)) + ~T^P( T ) ■ (V( r ) • ^tt±(^^))) + 



47tc|p(t)| c ' c 2m V 9 77 

+ f ^fliWr.flr)) + ■ („» ■ Ajf^r,^))), (B15) 



where 



A(r) = ^A ± (r,ff(r))- 
cm 

,2 



/Qmi i+i g Am Q m - 2/zm Q Amm,\ - 

" + H 7 2m" + 7 + 7 -^T) p{T) x 5(r ' r]{T)) - 

(B16) 



By evaluating the new momenta 



g ; (T) = g£lW = * + ^». (bit) 
^ .l^(mnf 



dr 

we obtain after some manipulations the Hamiltonian 



V / v dfj^r) , dff 2 (T) 
= \Jmjc 2 + (^(r) - Mr) ) 2 + + (^(r) - A 2 (r) ) 2 + 

+ 7^)1 " 7 p1r) - 7f±(T ' ^ (r)) " 7 ^7 * r) ' </ (r) ' ^) ^ (r ' ^ (r)) " 

Q Am Q m 2 - 2/im / 9 \ _ 

~ 7 — ^W(^)) " 7 2m 2 P(T) ' ( P(T) ' dif) 7r±iTMT)) - 

(B18) 

This can be brought to desired form using relative variables by substituting for and 
expanding the vector potential from the square root by using the Grassmann property. 
Consider first the non-relativistic limit (in which we keep terms of order 1/c). First we 

use 

yjm* c 2 + (jfrr) - Mr) ) 2 = M c 2 + nP(r) - -7^ == , • (B19) 

v y/mf c 2 + Kj 2 (r) 
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Note that we cannot use the Grassmann argument to eliminate the A (r) terms since 
they involve the nonzero Q 2 = —2 Qi Q 2 terms even with Q terms. We shall instead 
ignore them as being too high an order of 1/c 2 . So we get 



«i (r) • A (r) / « 2 (r) • A 2 (r) 



V 1 vX^ + ^tT) V 2 v ; v /m 2 c 2 + /5 2 ' 2 (r) 

+ 4^1^ " 7 p(r) ■ ^ (r ' ^ (r)) " 7 ^ p(r) " ( p(r) ' ^) " ±(r ' ^ (r)) " 

" I IT p(r) ^ (r ' ^ (r)) " I !!!! ^ p(r) ' ( p(r) ' ^) ^ (r ' ^ (r)) - (B20) 

In terms of the relative and collective variables 



iii = Wmf c 2 + ( — /c (r) + tt (t)) 2 m + 



m 2 c 2 + (^/?'(T)+7f'(r)) 2 



, , m 2 ^ ^ (^f/t'(r) -/(r)) -A(r) 

+ A /m 2 c 2 + ( — re (r) - 7T (r)) 2 m = + 



+ iMk " c p(r>7f±(r ' ^ (r)) - 7 ^ p(r) ■ ( p(r) ■ 4) Mt > ^ (r)) ~ 

Next, look at the limit as c — > oo in which 



L? C 2 " ^ - (r) + - (r)) 2 ^ m . c + 

y 777- 



2 



2mj c 

1 1 (^'(r) + R*+V(r)) 2 



v /m 2 c 2 + (^/5'(^) + (-) m ^(r)P 
and we keep only the terms of 0(l/c) and lower obtaining 



TOj c 2 to 3 c 3 



, (B22) 



+ - g ^ - g ^ ^ • fa) ■ §i) ^ ^ - 

- p(T)-TT ± (T,ff(T)) 2m 2 2 /im p( T ) • fa) • ^) *±(T,ii(T)). 

(B23) 

Now using Eq.(B16) we obtain in the non-relativistic limit 
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Hi c = mc 2 + 



+ 



K T 



m 



| 7f' 2 (r) 
2m 2/i 
g Q Am 



+ 



Q 



p(r) x B(T,ff(T)) A ± (r,r/(r)) 

c c m / c 



+ 



(-- + - „ S )n(r)-p(r)xB(T,ff(T)) + 



+ 



c 2m /i c 2m 2 fx 

® l( ^ 2 Q p(r)-7r ± (r,ff(T)) - Q ^ p(r) ■ (p(r) ■ -^j tt ± (t, ff(r)) 



4tt \p(r) 

Am _ m l — 2\im 

^rT^ ^ ~~ 2m? 



p(r)- (flr) — )Mt,v(t))- (B24) 



If we use Q^O, e = Q ~ <3i ~ — Q2> we get 



2 «V) e 2 

+ ^ + ^ • - P(r) x J3(r, ^(r)) + - ^ tt'(t) ■ p(r) x £(r, iffr)) - 
l\x m c c 2mfi 

- ep(r)-7f ± (r,?7(r)) ~ e ^-P( r ) ' (/( r ) ' 7^ J ^( r > ^( r ))> 

(B25) 

which agrees with Eq.(L.14) and (14.37) of Ref.[8] with e 2 = (one could also find the e 2 
terms of Eq. (7.2)). 



3. A Lagrangian for the Internal Motion of the Particles plus the Electro- Magnetic 
Field 

Let us consider the isolated system "charged particles with mutual Coulomb interaction 
plus the transverse electromagnetic field" inside the instantaneous Wigner 3-spaces. Its 
motion is determined by the Hamiltonian Mc = - S^ nt ) given in Eqs.(4.1) and then restricted 

— * — * 

by the rest-frame conditions V(i n t) ~ and by their gauge fixing }C^ nt ) ~ 0. Therefore we 
can look for the Lagrangian corresponding to the Hamiltonian Mc. 

By using the first half of the Hamilton equations with Hamiltonian Mc 
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dr 

dA±(r,a) o_ 
d7 



{fJi(T),M C } 



{A ± (T,a),Mc} = -7f ± (r,a), 



rrii c 



dr 



— ( df jj \ 

{dr) 



+ ^A ± (r,ff t (r)), 



2 C 



l-(^) 2 ' 



the inverse Legendre transformation leads to the following Lagrangian 



drjiir) 
dr 



dA 



7T± 



dr 



r,a - 



(B26) 



_ Wint) — 
C 



i 

Q1Q2 



87rc |?7i(t) - 77 2 (r)| c dr 



H ; A ± (t, a) 



dr 
+ 



+ 



1 

2c 



r dA A 
" dr 



= J d 3 a(22 s3 (^-'ni(r))[-m i c^ 



Q1Q2 



c \ or dr do 

dS 2 (r) 



dr ' 8nc \t}i(t) - rf 2 (r) 
+ 



dS 2 (r) 



dr 



dr 



(B27) 



In the last lines we have redefined the Lagrangian by isolating a total time derivative. 
This leads to the identification of the following generating function 
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S 2 = - C Y1 j d * a5 ^ B -Vi(r))Qiffi(T) ■ A ± (r,ff) 

= -y2Qifh(r)-A±(T,ffi(T)) = 



Q r 

— Pi2(r) • A ± (r,7fi 2 (r) + 
C L 



+ 
+ 

+ 



m 2 — mi 



Pu(r) ■ [Mr, , 
9 



9 



^±(r,77i2(r)) + 



?7i2(r) • (pi 2 {t) 



A±(r,^ 2 (r)) 



Q 

c 



i=l 



+ 



<9?7i 2 

^^(r).A ± (r,^(r)) + 0((p 

Q 

c 

m 2 — mi 

d 



12 " 



+ 

d , 
df) 12 ' 



A 



>Q*0 — \Mt) ■ A ± (t,T] 12 (t)) + 



m 



Mr) ■ (pi 2 (t) • ^-J A ± (T,ff 12 (T)) + 



+ Vu(t) ■ [ p 12 (r) ■ 



dff 12 



A±(r,ff 12 (r)) +0 (pia 







dff 12 



■d(r). A ± (r,ff 12 (r)) + 



m 2 -m 1 ( _ 



m 



d 



drfn 



+ ^Vn(r) • (p 12 (r) ■ -^-) A ± (r,ff 12 (r)) + o((p 12 • 



d_ 

df) i2 . 



A4r,Vn(r))\ + 

)• (B28) 



9 ^2 T 



drfn 



In Eq.(B28) we have shown the dipole approximation of S 2 . It contains an electric dipole 
term similar to the function S of Eqs.(B3) - (7.3) but with m-2 2 ^" 1 replaced by m2 ~ mi . 
Moreover it contains an extra term connected to the motion of the collective variable t] 12 {t) 
of the 2-particle system. 
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